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Abstract 

A new class of relativistic diffusions encompassing all the previously studied examples has 
recently been introduced in the article [l] of C. ChevaHer and F. Debbasch, both in a heuristic 
and analytic way. A pathwise approach of these processes is proposed here, in the general 
framework of Lorentzian geometry. In considering the dynamics of the random motion 
in strongly causal spacetimes, we are able to give a simple definition of the one-particle 
distribution function associated with each process of the class and prove its fundamental 
property. This result not only provides a dynamical justification of the analytical approach 
developped up to now (enabling us to recover many of the results obtained so far), but 
it provides a new general 77-theorem. It also sheds some light on the importance of the 
large scale structure of the manifold in the asymptotic behaviour of the Franchi-Le Jan 
process. This pathwise approach is also the source of many interesting questions that have 
no analytical counterparts. 
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1 Introduction 



The present article is at the confluence of two different stories that have met recently. 

The first was initiated by Dudley in a paper [2], written in 1966, where he describes the class of 
random Markov timelike paths in Minkowski spacetime whose laws are defined independently of 
any rest frame. These random paths represent the trajectories of particles whose speed is less than 
the speed of light, and whose laws are invariant by the action of the isometry group of the space. 
He proves in this article that there exists essentially a unique way of constructing random 
paths having the above properties. The phase space X H is well adapted to describe it. 

We write here M for the half-unit sphere {C = (t, x) € M x ; q{C) := - \x\l^^^ = 1, t > O} of 
the spacetime M x M^, equipped with the quadratic form q. The restriction of q to any tangent 
hyperplane of IHI is definite-negative. Any timelike path whose t-co-ordinate increases can 
always be re-parametrized in such a way that its speed belongs to M. Random timelike paths 
{7s}s^o ={70 + Jq ir dr^ are determined by their H-valued speed process {7s}s^o which has no 
other choice than being a Brownian motion on H (up to a constant time scaling). Minkowski 
spacetime thus has a canonical diffusion, in the same way as Brownian motion is canonically 
associated to Euclidean space. 

This fundamental work had to wait for the development of stochastic analysis and the article 
[3] of Pranchi and Le Jan, in 2005, to see its scope extended to the realm of general relativity. 
They defined a diffusion in any Lorentzian manifold using a stochastic development procedure 
similar in spirit to the construction of Brownian motion promoted by Malliavin and Elworthy, 
using stochastic differential equations in the orthonormal frame bundle of the manifold. 

The other story was born immediately after Einstein's theory of relativity and gravitation 
was accepted and spread in the scientific community. It deals with the extension of Boltzmann 
theory of gases to the relativistic framework. Although Boltzmann model is primarily a particle 
model of gases, most of the works have been on understanding the macroscopic behaviour of 
relativistic gases through the study of the raltivistic Boltzmann equation. One had to wait the 
nineties and the article |^ of F. Debbasch, K. Mallick and J. P. Rivet to see the introduction 
of a probabilistic mesoscopic model of diffusion of a particle in a fluid, under the form of a 
special relativistic counterpart of Ornstein-Uhlenbeck process. Generalisations of this model to 
the framework of general relativity have been given in later articles. 

These two stories have recently met with the proposition, made in the article [1] of C. Cheva- 
lier and F. Debbasch, to define a class of random processes including Dudley's process and the 
relativistic Ornstein-Uhlenbeck process, and characterized by the following property. There ex- 
ists at each (proper) time (of the moving particle) a (local) rest frame where the acceleration of 
the particle is Brownian in any spacelike direction of the frame, when computed using the time 
of the rest frame. The processes of this class were named relativistic diffusions in reference to 
the diffusion phenomenon they modelize. The authors of the article have started the study of 
this class developing an analytical approach to the situation based on a transport equation. We 
would like to propose in the present article a pathwise approach to this class of processes on a 
general Lorentzian manifold. With in mind the diffusion phenomenon of colloidal particles in 
fluids, we shall describe their dynamics as random perturbations of differential equations. In the 
spirit of the work of Franchi and Le Jan, we shall lift these dynamics to the frame bundle of the 
manifold, where they will be defined as flows of stochastic differential equations. This framework 
will enable us to re-prove directly many of the results obtained so far as well as new results and 
prospects stemming from the pathwise nature of our approach. 

We have organized the exposition as follows. Section [2] is dedicated to describing the class 



of relativistic processes in Minkowski spacetime, so as to separate probability and geometry 
problems. The class of relativistic diffusions is thus motivated and defined in section 12.11 We 
give in section 12.21 a probabilistic definition of the one-particle distribution function for each 
relativistic diffusion, and prove that it satisfies a fundamental equation. Section [3] is dedicated to 
investigating the general situation where the geometric background is any Lorentzian manifold. 
After having defined the dynamics in the orthonormal frame bundle in section 13.11 we shall 
spend some time in section [3^ looking at what can happen in the unit sub-bundle of the tangent 
bundle. We shall define in section [3.3. II the one-particle distribution function for each relativistic 
diffusion under a mild hypothesis on the global geometry of spacetime. The relevance of this 
notion in the study of the Poisson and Martin boundaries of the Franchi-Le Jan process will be 
discussed in section [3.3.21 Finally, we shall prove in section [33] a general i?-theorem. A number 
of open problems are scattered throughout the text. Numerous examples have been included so 
as to help the reader to get an idea of the state of the field. 

Notation. We shall write od for the Stratonovich differential. The sign d will be used for the 
usual differentiation with respect to the time, or for Ito's differential. 

2 Relativistic diffusions in Minkowski spacetime 
2.1 Definitions and examples 

a) Geometric framework. Recall Minkowski space is the product M x equipped with the 
metric 

VC = (i,x) € X M3, q{C) =e- ({x^Y + {x^f + (x^)') , 

if we write {t ^) for the co-ordinates of ( in the canonical basis {e°, e^, e^, e^} of M x M^. 

To distinguish Minkowski spacetime from the Euclidean space M^, we shall denote the former by 
E}'^. The half- unit sphere 

M = {C = {t, x) G R^'^ ; q{C) = 1, t > 0} 

inherits from the ambient (non-definite positive) metric q a Riemannian metric of constant cur- 
vature, which makes it a model of the (3-dimensional) hyperbolic space. As any timelike 
path can be re-parametrized so that its speed should belong to H, we shall look at the space 
E}'^ X H as the configuration space of timelike trajectories of a point of M}'^. The set of direct 
linear isometrics of q is the group 50(1,3). Any element g of 50(1,3) represents a rest frame 
g = (g°,g^,g^,g^) of M^'^. The function C € M^'^ ^ 9(g°,C) will be called the time function 
associated with the frame g. 

It will also be fruitful to define the motion of a(n infinitesimally small) rigid object. The 
configuration space of this dynamics will be the set M^'^ x 50(1,3). We shall look at a point 
(m, (g°, g^, g^, g"^)) as the infinitesimal rigid objeci0 m -|- ConvHull(5g^, 5g^, (5g^) contained in 
the affine spacelike hyperplane m -\- span(g^, g^, g'^) , and having 4- velocity g*^. An element of 
M^'^ X 50(1,3) can also be seen as an observer. 

Notice that 50(1, 3) has 4 connected components; we shall denote by 50o(l, 3) the connected 
component of the identity. To shorten notations, we shall write OM^'^ for M^'^ x 50o(l,3). 

The introduction of the following notations will clarify the description of the dynamics we 
are interested in. We shall denote by Ei £ so(l,3) the Lie element such that exp{tEi) is the 



^5 is some infinitesimal positive number. 



hyperbolic rotation of angle t in the 2- dimensional plane generated by e° and e*. In matrix 
notations 
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Four vector fields on 



^'■^ will be of particular interest. 

i/o((m,g)) = (g°,0), 

forz = 1..3, y,((m,g)) = (Cg^O- 



(2.1) 



Note that the M ' -part of the integral lines of the vector field i^o are the geodesies of . 
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{(m,C) G K^'^ X M} and write 
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which are straight lines. We shall set II 
{(m,g) eORl'3;gG50o(l,3)}. 

An important feature of our approach to relativistic diffusions is that we have chosen to 
describe the dynamics in the phase space OM^'^, where it has a natural and simple form; this 
corresponds to look at the motion of a small rigid object. We shall look at what happens in 
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in a later section. 



b) Dynamics. • Unperturbed system. We have indicated in the introduction that rela- 
tivistic diffusions should be considered as a class of toy models of diffusion in different media. 
We are going to define them as random perturbations of deterministic evolutions given by the 
flow of a vector fields V on M^'^. With in mind diffusion of particles in a fluid, we shall make the 
hypothesis that V has no M^'^-part and acts only on the 50(l,3)-part of OM^'^, although this 
assumption could be relaxed. The unperturbed sytsem is defined by the differential equation 

drris = ^„ds, 

' (2.2) 
dg.s = V{^s)ds. 

Note that the requirement that = g^ G EI implies that the parameter s is the proper time 
of the timelike path {ms}s^o of M^'^. 

• Action of the surrounding medium. How should we model the form taken by the 
random perturbation of the dynamics associated with a given medium? Maybe the proper way 
to proceed would consist in giving first a description of the microscopic thermodynamical and 
electro-magnetical properties of the medium in order to put forwards the source of randomness, 
and to infer from this description a description of the random perturbation it induces on the 
dynamics of a test object. We have chosen to propose a rather general action model which should 
convey the essential features of many situations, and not to model the medium itself. 

The action of the fluid on the moving object {e^js^o = { ("'■s, gs)}^>Q will be represented by 
the datum of an OM^'^-valued premsi6/e process {jsjs^o such that 3s(e.) = 3^ (("i.> g.)) = {"nis^fs) 
for some orthonormal basis fs = (/^ , , fg, ff) ofTm^M}'^^). The random perturbation induced 
by the medium on the dynamics results in adding to the deterministic acceleration a random 
part which is determined by the following requirement. When computed in the rest frame 3^, 
i. e. using its associated time, the acceleration of rUg has a deterministic part and a random part 

^Note that is and 63 have the same R^'^-part equal to nis. 



which is Brownian in any spacelike direction belonging to span(/^^, /^). To complete this 
description, we shall ask the vectors gl, , gs to be transported parallelly along the "Brownian" 
increment of g^. 

c) A preliminary example. Before giving a mathematically clean definition of this class of 
processes, let us look at the heuristic description of what happens when V = and the 'vertical' 
action process 3. is constant, equal to Id, i.e. fs = {e^, ••■,e^} for any s. 

Denote by {(m^, gs)}^^^ the OM^'^-valued process corresponding to these data and write 
for the e°-component of mg. As we have drus = g^ds, the function s i—> ts is a increasing 
function that can be used as a parameter of the process. Given t G M, set Tt = inf {s ^ ; tg = 
and look at the re-parametrized process {("in, gTj}t^g(^o^^Q); denote it by {{ffit,gt)}t^g(^^o^^^y 
The above description of the action of the surrounding medium on the dynamics means that the 
span(e^, e^, e^)-part of dg^ is a Brownian increment. 













— 

\ // 






// /T 



»1,3 



£1 



Figure 1: Dynamics when 3 = Id and V = 
The Brownian spacelike part e^odw\ of the increment of the speed can be seen in figure 1, 



j=1..3 



in red; the increment itself is in green. The notation w stands here for a 3-dimensional Brownian 
motion. If we write odg^ = g^od/J^, then 

i=1..3 

°dw\ = - q{e\gl)od(3i. 
j=1..3 

Denote by A{g) the 3x3 matrix with coefficients € [1,3]^ equal to q(e*,g-'). This matrix 
being invertible, 

odPt = -A{gty^odwt. (2.3) 
Back to the proper time s of the process, we shall write odg^ = g^od/?^. Write Ag for 

i=1..3 

A{gs). Identity (fOl implies that 

odps=q{e'',g^y2Aj^odws 

for some 3-dimensional Brownian motion w. The M^-valued process f3 is the process that really 
drives the dynamics. Last, we shall ask the vectors gjjgsjgs to be parallelly transported along 



the paths {g^ls^o in H- The above heuristic description gives rise to the following equations of 
motion 

odmg = g^ds, 
odgs = gsEiodPl. 

d) Definition. We shall now leave appart this example to write down the equations of the 
dynamics of {(m^, gs)}^-^g corresponding to general data V and 3. Recall the surrounding 
medium will be represented by the datum of a previsible process {is}s^o = {3s(e.)}sj>o such 
that 3s = (ms,/s) = (m^, (/^ , /g )) belongs to Omsl^^'"^- Its action on the dynamics has 
been heuristically described in paragraph b). Define the random matrix process {^s}s^o, with 
coefficient G [1,3]^ equal to q{fl,gi) at time s; set 

od(3s=q{flg's)'^Aj'odw,. (2.4) 

Definition 1. Define the -valued process (3 as above. A (V,^)- diffusion is a process 
{es}s^o = {("^«' Ss)}s>o stochastic differential equations 

odrus = g2 ds 

(2.5) 

odgs = V{es)ds + gs Ei odPl, 

where Einstein's summ,ation convention is used, as in the sequel. 

Using notations (|2.ip . equation (|2.5p can be written 

ode, = Hoie,)ds + y(e,)ds + V-(e,) odp^. (2.6) 

In reference to the interpretation of OM^'^ in terms of infinitesimal rigid objects given in para- 
graph a), this equation can be interpreted as describing the random motion of an infinitesimal 
rigid object in M}'^; there are nonetheless no need to understand it that way if you do not feel 
comfortable with infinitesimals. In any case, the simple and intrinsic character of this equation 
should be compared with the co-ordinate approach proposed up to now, as presented for instance 
in the article [5] of C. Chevalier and F. Debbasch. The simplicity of the formalism of stochastic 
differential equations will enable us not to rely on the covariant treatment used so far. 

Note that since 3s (e ) might depend on the whole history of e, until time s, the increment 
odPs shares this property, and equations (12. 4p and (12. 6p do not generally define a Markov process. 
This might be relevant from a modelization point of view if we consider an object having internal 
parameters evolving with time, and whose value at proper time s could influence the way the 
surrounding medium acts on it. Let us give three (Markovian) examples before commenting any 
further. 



e) Previously studied examples. Three (F, 3)-diffusions have attracted attention up to now. 

1. The Dudley (-Pranchi-Le Jan) process introduced by Dudley in [2] (and generalized in [3] 
by Franchi and Le Jan) is a perturbation of the geodesic flow. It corresponds to taking 
V = and 35 = 6,. The dynamics driving process /3 is then equal to the Brownian motion 
w, and no time-change is needecH. It is described in a simple way saying that 

• {Ssls^o is a Brownian motion on the hyperbolic space H, 
^That is, the time scaling q{fs, ga) is here equal to 1. 



• (gl,gs,gs) S obtained from (go, go, go) by parallel transport along the path 
{grlo^r^s, and 

• nig = mo + Jq Sr dr. 

This process is the only process determined entirely by the datum of the geometric back- 
ground (a result due to Dudley in [2](|fl))- This property gives it a special position in 
the family of (V, 3)-diffusions. Yet, its drawback as a model in Minkowski spacetime of 
a diffusing particle is that, except if we locate the source of motion in the particle itself, 
it is not clear what entity could give rise to such an interaction process. So it might be 
less interesting from a modelization point of view. Consult yet the article [6] of Dowker, 
Henson and Sorkin for a physical motivation from quantum mechanics. Nevertheless, the 
long-time behaviour of this process and its Lorentzian version may have many things to say 
about the geometry at infinity of spacetime; this might happen to be of some (theoretical) 
physical interest. We shall discuss this point in section [3.3.21 

2. The relativistic Ornstein-Uhlenbeck process (R.O.U.P.) was introduced by F. Debbasch, 
K. Mallick and J. P. Rivet in the article |^ as a model of diffusing particle in a fluid at 
equilibrium. It corresponds to the (y, Id) -diffusion with 



for some positive constant a. We have written here 7 for q{e^,g^) and grad for the gradient 
in H. In this case, the dynamics driving process f3 is not equal to the Brownian motion w. 
The existence for this process of an invariant measure of the fornil ae~^dm (g) dg found 
by Jiittner in [7] was a motivation for its introduction; see the introduction of the article 
[4]. We shall see in the general framework of section [321 that this OM^'^-valued diffusion 
gives rise to an HM^'^-valued diffusion, which is the R.O.U.P. as defined in [4] and the 
subsequent works of the authors and their co-authors. 

3. Last, Dunkel and Hanggi introduced in their article [8] a kind of mixing of the previous 
two models in which the frame 3^ = 65, as in the Dudley- Franchi-Le Jan diffusion, and V 
is constructed in such a way that the process admits the same awaited invariant measure 
as the R.O.U.P. 

We shall come back to these models in the general framework of section [3l 

f) Non-isotropic medium. This way of defining (V, 3)-diffusions has the advantage to be 
flexible enough to provide models of what should be a relativistic diffusion in a non-isotropic 
medium. We shall take into account the non-isotropy of the motion replacing the up to now 
isotropic input odws by a non-isotropic semimartingale in equations (|2.4p and (12. 5|) of dynamics. 
Setting for instance M = diag(l, 1, 2) and denoting by {Bs}s^o an M^-valued Brownian motion, 
the use in the R.O.U.P. dynamics of an input odvog = ModBg will give rise to a motion in a 
medium where one spacelike (fixed) direction differs from the others. One could also replace w 
by any continuous semi-martingale to adapt the model to a given situation. Jumps could even 
be introduced to take into account possible shocks. 

The article [9] of J. Franchi and J. Angst proposes another model in Minkowski spacetime of 
random dynamics in a non-isotropic medium. 

■^Note that we have uniqueness up to a time scaling by a constant in the H-Brownian motion {gsjs^o- 
^The measure dg is a Haar measure on the unimodular group S'Oo(l,3), and a and b are positive constants. 
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g) Probabilistic matters. Let us be more precise in the definition of a (y,3)-difrusioiil. Let 
{W, {'Ht}t^o) denote the Polish space C(M+, OM^'^) , endowed with the filtration generated by 
its co-ordinate process. Let 3 : 1R+ x W —^ OM^'^ be a previsible path functional. A {V,})- 
diffusion will consist in the datum of a filtered probability space (fl, {J^t}t^o,^) satisfying the 
usual conditions, an ({jTiji^O; IP)-Brownian motion w on M^ and a C(M+,OM^'3) -valued process 
e defined on (Q, {J^t}t^o) such that equations p. 41) and (|2.5p hold. These sorts of details will be 
implicit in the sequel. 

Existence and uniqueness results exist for equations such as (|2.4p and (|2.5|) . Consult [10] and 
the references given therein for example. These issues will raise no problem in the example we 
shall consider. 

We should apologize for the mis-use of the word "diffusion" in this context, as it is usually 
used when 3s(e.) = 3(6^), which is not supposed here. We have chosen to keep this denomination 
in reference to the situation it modelizes. The word "diffusion" will keep in the sequel its usual 
meaning, and we shall always write (y, 3) -diffusion for a process of our class. 

Last, we shall use the notation {es}s^o, indexing the trajectories by M+, regardless of the 
possibly finite lifetime of the process. One can add a cemetery point to the space to deal with 
such issues. 

2.2 One-particle distribution function of Markovian (F, 3)-diffusions 

As explained in the introduction, the main aim of his article is to convince the reader of the 
usefulness of a pathwise approach to relativistic diffusions. This section will illustrate this point 
giving a clear definition of the one-particle distribution function of a (V, 3)-diffusion. We refer to 
the article [H] of F. Debbasch, J. P. Rivet and W.A. van Leeuwen for a physical discussion of this 
concept of statistical physics and for the interest of a clear definition of this notior0. We shall 
investigate the general situation on a Lorentzian manifold in section 13.31 Let us first describe 
the framework of the problem. 

a) Framework. We shall suppose in this section that 

3s(e.) =d{es) 

for some function 3 : OM^'^ ^ OM^'^ such that 3((m,g)) = (m, (/°(e), /^(e))^ It follows 
that the process {e^js^o = { ("^si gs)}^>o is an OM^'^-valued Markov process. Write ^(e), or 
simply A, for the 3x3 matrix with coefficient G [li3]^ equal to ^(/^(e), g-'). The generator 
of the (y, 3)-process is given by the formula 

L:=Ho + V + ^ViB'^Vj, (2.7) 

where B = is a 3 x 3 non-negative symmetric matrix. Here as in the sequel, a vector 

field is seen as a first order differential operator; so, an expression like ViB^^Vjf should be more 
properly written Vi(^B^^Vj{f)) . The use of the notation 

A:=g(/°(e),gO) 

will be useful to shorten formulas, here as in the sequel. Recall we have supposed that the fiow 
of V preserves each fiber of the projection (m,g) € OM^'^ — > m G M^'^. 

^Refer to the chapter V.8 of the book |10) by Rogers and Williams for all this paragraph. 
'^The article [12] of W. Israel can also be consulted on this subject. 



We shall denote by dg the Haar measure on (the unimodular group) S0o{l,3) whose image 
by the projection g G 500(1, 3) i— > g'' is the Riemannian measure on H. Last, we shall associate 
to any subset A of E}''^ the (principal) bundle 

OA := {(m',g') € m' £ A, g' £ 50o(l,3)}. 

If A is a spacelike hypersurface of M^'^, denote by aA{dm') the volume measure induced by q on 
A; we define the measure 

YoLoA{dg' A dm') := dg' ® aAidm') 

on the bundle OA. 

b) One-particle distribution function. A few more notations are needed to define the one- 
particle distribution function and state its main properties. Fix a point e = (m, g) G OM^''^, and 
define the hyperplane of M}'^ 

Ve = {m e M^'^ ■,m'Gm + (g°)"^}; 
denote by the hitting time 

i?e = inf{s > 0; e, G OVe}. 

We shall associate to any a G 50o(l,3) and any t G M the hyperplane Yt^a '■= {m' G 
M^'^ ; g(m', a°) = t} and the hitting time Ht,a = inf{s > ; G OYt,a}- 

Notice that the Liouville measure induced by q on 0]Ri'3 is the product measure 

VOL((ig A din) := dg (8) LEB4(dm). 
We shall denote by L* the L^(VoL)-dual of the operator L; we have V* = —Vi and Hq = —Hq. 

Theorem/Definition 2. 1. Let Gq G OM^'^ he different from e. The random variable 
^He^He<oo has under Pep a smooth density /(eo; {m',g')) with respect to the measure 
VoLoVe(^^g' A dm') on OVg. 

The function e G OM^'^\{eo} i-^ /(eo ; e) is called the one-particle distribution function of 

the (y, 3)-difrusion started from eo. 

2. We have 

Eeo[/(eH,j] - J /(e)g(«",g")/(eo; e)VoLov,,„(de) (2.8) 
for any bounded function f on OYt,a- 

3. The function /(eo ; •) satisfies the equation 

L*/(eo;-) = (2.9) 

on OM^'3\{eo}. 

It is clear from its definition that this function is defined in an intrinsic way; physicists use 
to say that /(eo ; e) is a Lorentz scalar. We shall prove in section [3.3.11 a similar theorem in the 
general framework presented in section [3l We have chosen to present here a heuristic proof of 
point 2 and to give a detailed proof of the general statement after proposition [U in section [3.3. 11 
Points 1 and 3 of theorem/definition [2] are proved in detail below. 




Figure 2: Re-parametrized process 

< Proof - 1. The strategy of the proof is simple. Given e = (jn, g), we are going to re- 
parametrize the process as a function of the time associated to the frame g and see that 
/(eo ; •) is the density with respect to VoLoVe of the position at some fixed time of a hy- 
poelhptic diffusion. 

Define the chronological past o/OVe as the set /^(OVe): 

I (7(0)5 g') ^ OM^'^ ; 7 future-oriented timehke path from 7(0) to a point of the set m + (g^) 

g'G50(l,3)}. 

The random variable en^'i-HeKoo being identically equal to if eo does not belong to the 
chronological past of OVg, we shall suppose in the sequel that eo belongs to it, in which case 
He is almost surely finite. 

Set = q{mo,g) and define the stopping times 

VtGM, 5f =inf{s ^ 0; g(g°,m,) =t}. 

The process {est}t^t^ the process e. re-parametrized by the time associated with g. It 
has generator 




where Ag = Ag(g') = q{g^,g'^). We shall write 

esf = ((tg° + xsf),gs^) e OM^'^ with xsf G span(g\ g^ g3), 
and shall look at 

ef := {xsf,g'sf) G Ospan(g\g2,g3). 

The random variable is equal to mo -|- hg^ + ef^, with ti = q{m — mo,g^). We shall 
prove the first point of theorem [2] showing that the Ospan(g^, g^, g^)- valued diffusion e^ is 
a hypoelliptic diffusion. The distribution at time ti of this diffusion will then have a smooth 
density with respect to the volume element on Ospan(g^, g^, g'^) to be defined below. The 
measure VoLoVe being the image of the volume element by M^'^-translation by mo -|- hg^, 
this will imply that en^ has a smooth density with respect to VoLoVe • 

To complete this program we shall denote by := (x, g') a generic element of O span(g'^, g^, g^ 
Note that since g has determinant equal to 1, the change of variable formula says us that 
the volume element induced by q on the 3-dimensional vector space spanned by g^ , g^ , and 



is the 3- dimensional Lebesgue measure. We shall write VoL(Q)gp(gi g2 g3)(dg' A dx) = 
dg' (8> LEB3((i2;) the volume measure on the bundle Ospan(g^,g^,g^). 

To describe the generator of the process e^, denote first by dx the differentiation operation 
in the direction of the vector space span(g^, g^, g^) and decompose g"' as 

g'°=Agg°+x:i.g^ 

i=l 

Using these notations, we can write for any smooth function / 

Ag Ag Ag /Ag 

So the generator of the process is given by the formula 

Ag Ag ^Ag 

Write heo{t,dsi^) for the law of ef,t ^ Iq. As is well known, these distributions satisfy the 
heat equation 

dthl = {L^y^h%, (2.10) 

where (L^)*^ is the L^(VoLfjgp(gi g2 g3-))-dual of the operator L^. Since the matrix B = 
(^A~^YA~^ is symmetric and V*^ = —Vi, we have 

{L^Y'hi = - {dxhf^) (^) + y-(^) + l{v.B^%) (a^). 

It is easy to see on this formula that the operator dt — [L^)*^ on MxO span(g^, g^, g^) satisfies 
Hormander's criterium for hypoellipticity. It follows that hfg{t, •) has a smooth density with 
respect to the measure VoL^gp^gi g2 g3)((iaS) on O span(g^, g^, g^), for t > g(mo,g'^). We 
have seen that it implies that e^^^ has a smooth density with respect to VoLoVe • 

2. As said above, we present here a heuristic proof of point 2. The reader will find the 
detailed proof of the general statement after proposition [H in section 13.3.11 We are going 
to explain the situation for Dudley's process, nothing else than additionnal notations being 
necessary to understand the general case of Markovian (y, 3)-diffusions. 

• We shall get a clearer image of the situation considering the continuous dynamics described 
by equation 

des = Hods + Viodwl (2.11) 

as the dynamics of a random walk {es}s^o = {{fhs,gs)}gyQ making infinitesimal steps. 
Given an 'infinite' integer N {i.e. a nonstandard hyperfinite integer), the quantity ;^ is a 
positive infinitesimal. Let us denote by {A^} a 'sequence' of iid M'^-valued centered Gaussian 
random variables with variance The dynamics of the random walk is defined on each 
interval of the form [jj, ^^), A; ^ 1 as follows. 

• The process gs has a jump at time j^: g± = g/j^^ - . exp(£^j Ait;^) . The process 

^ N i JV j 

{"T-sjs^o has no jumps at that time. 

• gs is constant and drhg = g'^ds, in the time interval (-^, ^^). 



gs is constant and drris = g^fis in the first interval [0,-^). The statement "The random 
walk {es}s^Q provides the solution of equation (|2.1ip " can be given a precise meaning in the 
framework of non-standard analysis, and holds true, when correctly interpreted. This way 
of saying things is, in any case, useful (justified) and intuitive. 

Notations. We shall denote by Peo the law of the random walk started from eo. Given two 
possibly infinite real numbers a and b, we shall say that a and b are equal up to a negligeable 
quantity if ^ is infinitesimally close to 1; we shall write a b. The notation Haar(-) will 
stand for a Haar measure on 5*00(1, 3). 

• Formula (|2.8p will hold true if we can prove it for any function / on ONt.a of the form 

/(m',g') = lA("/)lG(gO, 

for sufficiently small infinitesimal open sets A C Nt,a and G C 500(1,3). We shall suppose, 
without loss of generality , that ^ x G is a (connected) neighbourhood of a given point e = 
(m,g) € OVf^Q. We shall associate to e the hyperplane Vg = {m' € M^'"^ ; m' G m+ (g*^)"*"}. 
To distinguish the Lebesgue measures induced by q on q, and Ve, we shall denote them 
by LEBg" and Leb| respectively. 

If we let A' be the set of points of Vg of the form x + sg'", for a; G ^, s € M and g' G G, the 
LEB|-measure of A' is equal to 

LEBi(y')^'?("°,g°)LEB*'°(A). 

Let now M be an infinite integer and let run M independent infinitesimal random walks 
started from eo G OM^''^. Write NaxG and Na'xg for the (random) numbers of trajectories 
of the random walk that hit OVt^a and OVe in Ax G and A' x G respectively. If A is small 
enough for NLeb'^^"' (A) to be infinitesimally close to and M is large enouglo, (P^^j*^- almost 
surely) 'almost all' the trajectories of the random walks hitting A x G will hit it in a time 
interval where gs is constant. As the length of this time interval is much bigger than the time 
needed by any timelike path to go from A to A' , the trajectories of the random walk will hit 
^' X G on the same time interval where they hit A x G. As only a negligeable quantity of 
trajectories hitting A' x G will not hit A x G, we shall have on the one hand 

NaxG-Na'xG, Peo^'^ - almost surely. 

As the strong law of large numbers ensures us that 

NaxG X /i^^(e)Haar(G)LEB3'°(^) 

Na'xG X /(eo ; e) Haar(G)LEB^(^') ~ M x g(a°, g°)/(eo ; e) Haar(G)LEB3'"(^) 
on the other hand, it follows that 

<(e)^(?(a°,g°)/(eo;e). 
Both quantities being standard reals, we actually have equality. 

3. We are now going to use equation (12. 8j) to give a proof of equation (|2.9p . This will be 
done fixing a frame a G 50o(l,3) and proving that we have L*(^] = 0, where we have 
denoted by Aq the function e = (m, g) q{a^,g^). 

*Equal to an infinite integer depending on A'' and Leb3°'(A). 



A frame a having been chosen, define the stopping times 

Vr G M, 5,. = inf{s > 0; g(a°,ms) > r}, 

and the process e" in the same way as the process has been defined above. It can be 
proved as above that the random variable e" has a smooth density h^^^{r^ ■) with respect to 
VOL(Q)gp(„i ,^2 q3), under Peo! if is defined for r ^ qivtiQ^c^). The function h%^(r^ •) is defined 
as equal to for r < q{mQ, a^). Identifying (r, (m', g')) G M x (Osp(a^, a^, a^)) to the point 
(rQ° + m',g') of OM^'^, the function h'^^ will be seen as a function on OMi'3\{eo}. Three 
more notations will be needed: will stand for the generator of the O span(a^, a^, a^)- 
valued diffusion e", we shall write Dmh^g for the partial differential of h^^ with respect to 
m([f|) and use the notation to refer to the partial differentiation operation in the direction 
of span(a^, a^, a^); last we shall decompose a vector g*^ G BI as 

3 

gO = 5(gO,aV + ^i;,a\ 



Note the relation 



i=l 



{DmhZ) ( 1 = drK. + (D^K 



which can be written 

- {^^"1) (^) = (^) + (2.12) 

Recall that we write for q{a^,g^). It can be proceeded like in the proof of the proposi- 
tion/definition [2] to show that /ieo('i ') satisfies the heat equation 

drK, = {L^Y''hl, (2.13) 

where 

and the operation is the L^(VoL(Osp(„i^„2^„3))-dual operation. Using equation (I2.12p . the 
heat equation (|2.13p can be written 

-Ho[f^)+V*^[^)+-{V,B'^V,)[\f^)=0. (2.14) 

Note that since the vector field V acts only on 5*00(1,3) we have y*" = V*; we have 
recalled above that Hq = —Hq. So, equation (|2.14l) can take its final form: L*(^j = 0, 
i.e. L*/(eo; = 0. > 
This theorem/definition needs a few comments. 

• Formula (|2.9p is fundamental in the approach developped by Debbasch, Rivet and their co- 
workers. Their analysis of the situation entilery rests on a similar transport equation. Although 



^The map D^/i" , is for any e £ OR^'^ the linear form ( G K^'^ ^ iiin!!mllt!i£L_!!silll; this limit is denoted 
1^ 77,0 ^ 

by (^™/iSo|e)(0, or simply {D^KJiC). 



it can be argued that since equation (12.8(1 implies that the one-particle distribution function 
determines the hitting distributions of the process at any times of any rest frame, a theorem of 
Blumenthal-Geetor and McKean ensures us that this function essentially determines the process, 
such a position should be taken with care. Indeed, the development of stochastic analysis has 
shown that one can gain much insight in the situation looking at the pathwise behaviour of 
processes rather that looking at analytic quantities such like their semi-group. We hope to 
illustrate this point throughout this article. In any case, theorem [2] makes it clear that the 
fundamental quantity is not a hitting distribution h'^^ but the one-particle distribution function; 
a fact which was not put forwards in the article [I^ of C. Barbachoux, F.Debbasch and J. P. 
Rivet. 

• Equation ()2.9p has a clear meaning from a Markov process point of view. It says that the 
measures /(eo ; e)VoL((ie) on OM^'^ are invariant for the (F, 3)-diffusion. It is tempting to ask 
wether these measures and their possible renormalized limits as eg goes to infinify are sufficient 
to describe the set of all invariant measures. For instance, it would be interesting, in the study 
of the R.O.U.P. in Minkowski space, to see if the strong recurrence of the process {gsjs^o is 
sufficient to prove that the Jiittner measure ae~'^'\o'L{de), alluded to above, is the only measure 
we obtain sending eo at infinity, while imposing the limit measure to have mass in any open 
sei0. Even though a complete answer of the general question is out of reach at the moment, we 
shall come back in section 13.31 to related matters in the general framework that we are going to 
present now. 

3 Relativistic diffusions in a Lorentzian manifold 

We shall now proceed to defining (V, 3)-diffusions in a Lorentzian manifold. Let (M, g) denote 
a (1 -|- d)-dimensional Lorentzian manifold, endowed with its Levi-Civita connection. As in 
Minkowski space, we shall construct the dynamics in a bigger space than M. We shall first recall 
in section ISTTl a) how one can construct this space (the orthonormal frame bundle over (M, g)) 
and the analogue of the above vector fields Hq and Vi before defining the class of iV, 3)-diffusions 
in section [3]TJ b). We shall then see in section 1X2] that some situations give rise to a sub- 
diffusion in the (future-oriented) unit tangent bundle of M. Several example will be discussed 
before returning in section [3?3] to the study of (F, 3)-diffusions. We shall define in this section 
the one-particle distribution function of the (V, 3)-diffusion and prove its fundamental property. 
This result will shed some light on the structure of L-harmonic functions (section [3.3. 2p and will 
provide a simple proof of a general H-theorem (section [3^ . 

Hypothesis. We shall suppose from now on that (M, q) is oriented and time-oriented. 
3.1 (V, 3)-diffusions in OM 

a) Geometrical objects in play. Given some point m € M, it will be useful to consider an 
orthonormal basis {g*^, g'^} of the tangent space r^M to M at m as an isometry from (R^'^, q) 
to (TmM, so, strictly speaking, g* = g(e*). 

The orthonormal frame bundle of M is just the collection 

OM = {(m, g) ; m e M, g an orthonormal basis of TmM). 

^"Other measures can be obtained if we do not impose this condition. 
^^The letter q has here two different meanings. 



We shall write OU = {(m, g) ; m e Z//, g an orthonormal basis of TmM} for any subset U of M. 
One defines the manifold structure of OM as follows. This structure being local, it suffices to 
define the structure of OU for any (small) domain U of M; take it small enough to be the domain 
of a chart x : U ^ M}'^'^. Applying Gram-Schmidt orthonormahsation procedure to the family 
of vectors {S^^jj^Q 3 in each tangent plane, on defines a section cr : W — > OU. The identification 

i:U X 0{l,3) ^ OU, {m, g) ^ {m, a{m)g) 

gives OU its differentiable structure (compatible with changes of chartsjll. Note that 0(1, d) 
acts on OM on the right: the action of g' on (m, g) in the above chart i is 

{m,g).g' = {m,a{7n)gg'). (3.1) 

Note that OM has several connected components. We shall be interested in dynamics leaving 
these components globally fixed. We choose to consider only one of them, specified by the 
requirement that g'^ should be future-oriented and that the orientation of g should be direct (we 
have supposed the space oriented). The above action of the connected component of identity in 
5*0(1, d) preserves our connected component. We shall also denote it by OM, as there will be 
no risk of confusion. 

Action (|3.ip enables us to define vector fields on OM: 

Vi{{m,g)) = ^ ((m,g).e*^0, i = l..d. 
dt\t=o 

Last, we shall define the vector field Hq as the infinitesimal generator of the geodesic flow on 
OM. The dynamics {(m^jg^)} of this flow is described by asking that = g^, and gs should 
be transported parallely along the path {rus}- One has for instance i/o(("^,g)) = (g'^,0) in 
Minkowski's flat spacetime, in accordance with the previous definition of Hq given above. 

Notation. We shall write e for a generic element of OM. 

b) (F, 3)-difrusions. We are going to define (F, 3)-diffusions following the same approach as in 
Minkowski space. We shall thus consider these diffusions as random perturbations of the flow of 
a differential equation in OM of the form 

deg = Ho{es)ds + V{es)ds 

where V is any vector field on OM. As in section 12.11 we shall not modelize the surround- 
ing medium itself but just its action on the dynamics. This action will be given through the 
datum of an OM-valued previsible process {^s}s^o such that 3s(e.) = 3s((W'.,g.)) = irns,fs) 
for some orthonormal basis fs of Tm^M.. Roughly speaking, it has the property that, when 
computed in the rest frame 3s(e.), i.e. using its associated time, the acceleration of m. has a 
deterministic part and a random part which is Brownian in any spacelike direction belonging to 
span(/i(e,),/2(e,),/3(e,)). 

Define {As}s^o as the d x d random matrix process with coefficient (i,j) G [Ij*^]^ equal to 
QifhSs) at time s, and set 

odp, = qif^,g'S^A-'odWs. (3.2) 



Consult for instance chapter 10 of the book [14] of P. Malliavin. 



Definition 3. Define the W^-valued process 13 as above. A (V,^)- diffusion in (M, g) is an 
OM-valued process {e^j^^o = §s)}s>o ^^^y^fy'''''^9 stochastic differential equation 

odes = Ho{es)ds + V{es)ds + ^^(e,) odPl (3.3) 

If you do not feel comfortable with this stochastic differential equation, we shall give a step- 
by-step description of the dynamics in the next section. The remarks on probabilistic formalism 
and existence and uniqueness results made in section 12.11 g) apply here. Let us emphasize the 
interest that the above general definition might have for modelization. It provides a model of 
evolution of an object which has internal parameters (such as a spin) influencing the way it 
interacts with the surrounding medium, and whose value at some proper time depends on its 
past history. Challenging questions arise from this non-Markovianity of the model; yet, as only 
Markovian examples have been studied so far, we shall mainly explore this situation in the sequel. 

Example: Franchi-Le Jan diffusion using co-ordinates. This diffusion is the (0,e )- 
diffusion, first defined in [3j. Note that the (0, e.)-diffusion is essentially the unique {V, 3)-diffusion 
determined entirely by the geometric background (M, g). We asked in section \27\\ e) which en- 
tity could give rise to the random excitement Vi{es) odw\ that enters in the equations of motion 
of the Dudley(-Pranchi-Le Jan)-diffusion process in the empty spacetime of Minkowski. This 
objection disappears when we consider the (0, e.)-diffusion in any spacetime (M, g) containing 
matter. It is in that case possible to add to the macroscopic description of matter given through 
the stress-energy-momentum (non-null) tensor a microscopic (quantum) description of matter 
from which randomness can be infered to comJ^. 

Equation (13. Sp takes for this process the form 

odes = HQ{es)ds + Vi{es) odw\. 



To describe how we can write equation (|3.3p using co-ordinates, note first that the data of local 
co-ordinates on M induces local co-ordinates on TM: a vector p G T^M will be written 
p = ^ V^dxi- Denoting then by T : ]R^+°' x M^+'^ ]R^+'^ the Christofel map associated with 

1=0. .d 

these co-ordinates, the dynamics of the Franchi-Le Jan diffusion takes the form 

odvris = ggds, 

odgO = -r(gO,gO)(is + glodw's, (3 4) 

j=0..d 

odgi = -r(g°, g^) ds + g°s°dwi, for j = l..d. 

These equations have to be written using the preceding co-ordinates. If one wishes to use Ito 
differentials, the system becomes 

dnis = g°ds, 

dg°= f-r(g°,g°) + ^g°l ds+y^ gldw] 



2 v ' (3-5) 

d^s = (-r(g°, gD + ds + g° dwi, for j = l..d. 



^■'Consult the article [6] for results in this direction. 



Remark that if we write Qm the matrix of the metric in these co-ordinates at point m, then the 
co-variance matrix of the martingale g* dwl is equal to gs(gs)* — The fact that it 

i=0..d 

depends only on nis and implies that the sub-process {(ms,gs)}g>Q is itself a diffusion. The 
investigation of such situations is the object of the next section. 

Note, en passant, that since we can read the matrix Qm^ on the co- variance of the martingale 
part of gg , it means that the local geometry of (M, q) can be recovered from the pathwise study 
of the sub-process {(m^, g^)}^^^. To determine what amount of information on the large scale 
structure of the space (M, q) one can obtain from the pathwise study of this process or of the 
(y, 3)-diffusion is a much harder task; we shall come back to it in section [3.3.2[ 

The heuristic explained in section [2?T] and motivating the above definition of (F, 3)-diffusions 
should make it clear that (F, 3)-processes should be considered as models of diffusion in a ho- 
mogeneous medium. Note yet that the input of a non-isotropic excitement in place of dw in 
equations (|3.2p . (13.31) would provide models of diffusions in a non-isotropic medium. 

Hypothesis for the remainder of the article. With in mind the diffusion of particles in a 
fluid, we shall suppose from now on that the flow of the vector field V leaves each fiber of the 
projection (m, g) G OM i-^ m G M stable. 

3.2 Sub-diffusions in HM. 

As emphasized in section 12.11 a) in the framework of Minkowski spacetime, (V, 3)-diffusions 
defined above can be considered as models of random motion of an infinitesimal rigid object in 
a relativistic medium. It might be interesting in some situations to define what could be the 
random motion of a point in such a medium. To investigate a physically motivated classical 
framework, we shall concentrate on Markovian processes. 

As noted after definition [U the (y, 3)-diffusions are not Markovian unless we choose a Marko- 
vian previsible process 3: 

3s(e.) =3(es). 

This requirement is not sufficient yet to ensure that the sub-process {(ms,g^)}^^Q of {es}s^o 
is itself a Markov process. We give in paragraph a) a simple condition which is proved to be 
sufficient in paragraph b). Several examples are examined in paragraph c). 

Throughout this section, we shall suppose 3 regular enough to have existence and strong 
uniqueness in the system (|3.2I) , (13. 3p . We shall denote b\l^ 

MM = {(m, g°) € TM ; m € M, g° e TmM future-oriented unit vector} 

the (future-oriented) unit sub-bundle of TM. This space is the phase space of the set of 
timelike paths in (M, g). The map 

TT : OM M 

will denote the projection (m, g) 1-^ m, and vr : OM HM the projection {in, g) i—>- (m, g^). 

a) A sufficient condition to a have a sub-diffusion in HM. In addition to the hypothesis 
■K^V = made above, we shall suppose that 



We write here g'' for the vector of its co-ordinates in the basis {9^i}i=o..d. 
^'Recall that we have supposed (M, g) time-oriented. 



• there exists a vector field V on HM such that V is the horizontal lift of V to OM(0). 

We shall begin our investigation with the particular case of the (0, e,)-diffusion of Franchi and 
Le Jan. Remember equation (13.11) describing the action of 0(1,3) on OM. This action induces 
a right action of 0(3) C 0(1,3) on OM, which amounts to rotate the vectors g^,g^,g^ in the 
Euclidean space they generate and leaves HM C OM stable. Given the Brownian input w in 
equations (|3.2p . (|3.3|) . denote by 

e{s,eo;w) (3.6) 
the (unique strong) solution started from cq. We have for any ^£0(3) 

7r(e(s,eo5; w)) =7f(e(s,eo; gw)). 

Since gw is also a Brownian motion, the law of {7f(e(s,eo; 9'w))}gyQ does not depend on g G 
0(3), but only depends on tt (eg) € HM. The sub-process {tt (e(s, cq ; w)^ }s>o ^^^^ a diffusion 
in HM. 



b) A step-by-step description of the dynamics. The general case is covered by the fol- 
lowing theorem. 

Theorem 4. Suppose there exists a function 3° : OM MM such that 

• 3°(e) = {{rrijg)) G H^M depends only on (m,g'^), and 

. 3(e,) = (m„(30(e,),/^(e,),...,/(e,))) ^^'^^ 

for some functions f^, ...,/'^. Let {m,g^) € HM. Then, given any choice of G OM such that 
7r(eo) = (m, g*^), the law of the MM-valued process 7r(e(. ,eo; w)) depends only on {m,p) and 
the function j'^ , and not on the particular choice of f^ , and e^. The process 7r(^e[.,eQ ; to)) 
is a diffusion in HM. 

We shall present a heuristic proof of this fact, the remaining work being just a matter of 
formalism. Equation (|3.2p and (13.3(1 are the mathematical expression of the following heuristic 
dynamics explaining how on constructs Bg+Ss from e^. 

1. Set rus+ss = TTis + g° ds, 

2. then, set g^^.^^ = g° +(5g^ + V(m^ gO)(5s. The increment (5g° is the only vector of Tgo(HmM) 
such that its projection in span(/^(es), /'^(e^)) parallelly to i^dnisjg^)) is equal to the 

scaled Brownian increment '?(3°(es), g^) ^ y~^/^(es) odwl. 

i=l 

3. Last, transport parallelly {gj, •.•,gf} along the increment 6g^ + ^(m^.g^) of §«• 

Examining this description of the dynamics, we see that any previsible orthonormal transform 
of the basis {/^(e^), /"^(es)} will leave the law of the Brownian increment unchanged, so that 
the law of 6g^ will also be left unchanged. Note also that the changing eg G OM to another 
starting point with the same HM-projection will only influence the dynamics of g^, ...,gf . These 
remarks justify theorem HI To put this argument in a polished probabilistic form is a matter of 
formalism. 

"Denote by ipt{-) the flow of the vector field V on OM and by {'?t}t^o the flow of V on HM. The above 
hypothesis means that the point v't(m, (g", g^, g'')) G OmM is obtained by parallel transport of (g^,...,g'') 
along the path {'Ps{{m,g°))} in H^M. 



c) Examples. 

1. Dudley- Franchi-Le Jan diffusion in Minkowski spacetime ([2], [3]). We have 
already given its description in 12.11 f), 2: {g^j^^o is a Brownian motion on EI and = 
^0 + /o Sr dr. The usual stochastic development procedure can be applied to this process 
to construct its HM- version from its HM^'^-version; see [3], theorem 1. 

2. R.O.U.P. in Minkowski spacetime ([!]). This process is the HR-'^'^-sub-process of the 
(y, Id)-diffusion on OM, where y((m,g)) = — agrad(ln7), for some positive constant a, 
and 7 = q{e^,g^). In this flat spacetime with global co-ordinates {t,x), the dynamics 
may be re-parametrized by the time t; the state space then becomes {ix,q) eR^ X R^}, 
where {t,x) are the co-ordinates of m and q is the span(e^, e^, e^)-part of g". With these 



notations, 7 = 7(g) = yl + l^lly^j. The step-by-step description of the dynamics (or, 
more formally, the stochastic differential equation (|3.3p ) immediately yields the following 
stochastic differential equations for where w is an M^-Brownian motion: 

dxt = -j^dt, 
dqt = —2a— — rdt + odwt] 
this is the original description of the R.O.U.P. up to some constants. 

Notice that the process {qt}t^o is a Kolmogorov diffusion in R^. It has a unique invariant 
measure /i, which is a probability and has density with respect to Lebesgue measure pro- 
portional to e~^'^"'^'^\ As we have lim |g| ,^(|4a V7P — 4aA7)(g) = 16a^ > 0, a well known 
theorem ensures us that satisfies a Poincare inequality. As is also well knowrl^. this 
implies that the semi-group of the porcess {qt}t^o converges to equilibrium exponentially 
fast in at least like e"^^*-" -S)t^ S > 0. This fact sheds some light on the 

numerical simulations made in section 4 of the article [4]. 

R.O.U.P. in an arbitrary inertial frame0. It might be enlightening to write down 
the equation of the dynamics using the time r and the (x,q)-co-ordinates associated with 
any orthonormal frame g of R^''^. We shall write F(r, q^) for the damping force in these 
co-ordinates; note that is depends on r and p^. Its precise expression is unimportant. 

To take advantage of the description of the OM^''^ process given in the above step by 
step description of the dynamics, and to take advantage of the irrelevance of the precise 
orthonormal frame {gs,gs,gs} of TgoH we use in this construction, we chose to take as a 
basis of TgoM the family 

{g^ - (gSg°)g°, g^ - (g^g°)g°, g' - (g^ g°)g°} 

and write down the Vect(g^, g^, g^)-part of the increment of dg^ as 

(*) = E IE (g' - (g'' Ms')odpn gK 

fc=1..3 \i=1..3 / 

In this expression, the matrix Ag used to define (3 has coefficient (i, j) equal to (e*, g-' — (g-' , gg 
it depends only on g^. Switching from the description in terms of proper time s to the 

^'^See for instance the book |15) of C. Ane et al. 
^^Compare with the article |13| . 



1 , s ^1 

2 



evolution in terms of time r results in multiplying (*) by y-^j = q{e j gs) ^ = li^ir 
This finally gives 

dxr = '^^ dr, 

7(qr) 

dq, = F(r,q,)dr +7(q,)-i ^ ( ^ (g^' - (g^ gO)gO, g'^) 

A:=1..3 i=1..3 



(3.9) 



The vector g^ G H is determined by q,. No other choice of g^, g^, gf would give something 
fundamentally simpler. This complicated expression of the dynamics means nothing else 
than the inadequacy of the choice of co-ordinates to describe it. 

3. R.O.U.P. in the spacially flat Robertson- Walker spacetime (|16|). This model on 
expanding universe is the product R x equipped with a metric of the form dt'^ — a{t)'^dx^ , 
where a > 0. As in the preceding example, one can describe the trajectories of the R.O.U.P. 
using the absolute time t and the state space {{x,q) G M'^ x M"^}. We shall write 7t(g) = 



l + a{tf\q\ The step- by-step scheme (or equation (|3.3p) yields the equations of 
dynamics: 



dxf = — T-^dt 



it{qt) 

2 

ItiqtT" ' a{t) 



qt 1 ^^-^^^ 

2a a[t)^ — - — -dt -\ — — odwt; 



the gradient part in gives rise to the term —2a a(t)^ — '^^—dt, w is an M'^-Brownian mo- 

Itiqt) 

tion, and the —-^ odwt term is the Brownian increment in the Euclidean space (M'^, —a(t)'^dx'^^ 
a{t) 

It can be proved that this diffusion has an infinite lifetime. 

4. Franchi-Le Jan diffusion in the spacially flat Robertson- Walker spacetime. We 

shall use the notation (m, g*^) for a point of HM to describe the dynamics of this process. 
Using the canonical co-ordinates {t,x) in M x and denoting by [{t,x), {i,x)) the associ- 
ated co-ordinates in TM, we have seen in equation (13. Sh that the equations of the dynamics 
take the fornrl 



dnis = ds, 

dis = (^is - (aa')(t.)||x.|||„,i^ ds + dM*, 



dxs= [^-2^i,]±,ds + dM^, 



where the M'^-valued local martingale M = yM^,M^j has co- variance 

isXs XsX* a"2(ts)Id; 



^^Compare the derivation of these equations with the approach of the article |T6] of F. Debbasch. Note that 
the dynamics is described in this paper not in HM C TM but in T'M. 

^"Consult for instance proposition 35, p. 206 of the book |17) by O'Neill for the computation of the Christoffel 
symbols in a warped product. 



It has been shown by J. Angst that this diffusion has an infinite hfetime. Remark that 
the M^-valued sub-process {(ts,4)}g>Q is a diffusion. This kind of decomposition of the 
diffusion into smaller dimensional diffusions has been the key of the previous investigations 
in Schwarzschild and Godel's spacetimes. See [3] and the article [18] of J. Pranchi. 

3.3 One-particle distribution function 

The aim of this section is to clarify the so-called notion of one-article distribution function in the 
general framework of Markovian (F, 3)-diffusions on any Lorentzian manifold. We shall define it 
properly in section 13.3.11 and prove in theorem [TT] that it satisfies a remarkable equation. This 
theorem will justify the analytic approach developped by F. Debbasch and his co-workers, as 
exposed in [5] or [l] and the references cited therein. The relevance of this notion in the study 
of L-harmonic functions will be described in section 13.3.21 

The approach to one-particle distribution functions developped here is similar in spirit to the 
physical approach exposed in the article [H], in a physical/mathematical stylj^. It should be 
noted yet that only the special relativistic situation is investigated in this article, whereas we 
deal below with the general relativistic case. 

In order to ease the understanding of the situation, we shall make a hypothesis on the 
global geometry of the space (M, q). We shall suppose the spacetime(M, g) strongly causal: 
every point of M has arbitrary small (connected) neighbourhoods which no non-spacelike paths 
intersect more than once. This is a mild global assumption on the geometry of the space, satisfied 
by most of the models of physical spacetimes. This excludes, yet, pathological spaces where closed 
timelike paths exist, like Godel's spacetime. 

We shall also use the following local property, shared by all Lorentzian open manifold. Any 
point has an open (relatively compact connected) neighbourhood on which a time function is 
defined. By time function we mean a smooth function whose level sets are spacelike hypersur- 
face@. We shall denote by Um such a neighbourhood associated to a point m G M. We shall 
suppose Um small enough to have the property that no non-spacelike paths intersect it more 
than once. As a consequence, it will enjoy the following extra property. Any timelike path in 
(M, q) will hit any spacelike hypersurface of Um at most once. This property will be the main 
ingredient used to define of the one-particle distribution of (F, 3)-diffusions. 

3.3.1 One-particle distribution function. 

The initial point eo of the (F, 3)-diffusion will be fixed throughout this paragraph. Given a point 
e = (m,g) G OM, different from eg, define the collection 

Ve = {V; spacelike hypersurfaces of M contained in Um and such that m G V and r^V = (g^)"*"}. 

Associate to any V G Ve the hitting time 

H = inf{s ^ 0; jn^ G V}. 

Given any point m' in M we shall denote by VOLm'(^g) the Haar measure on Om'M, normalized 
in such a way that its projection on H^'M is the Riemannian volume element induced by q. 
Recall the definition of OV = {(m,g) G OM; m G V,g G T^M}. Let us insist on the fact that 
even if V is a sub-manifold of M, the element g' of a point (m, g) G OV is not an orthonormal 



Consult also the article |12) of W. Israel for a similar point of view. 
^^We shall construct these neighbourhoods in the beginning of the proof of proposition/definition O 



basis of Tj^V, but an orthonormal basis of T^M.. We shall write e = (m,g) for a generic element 
of OV and shall denote by a^{dfh) the volume element induced by q on V. With these notations, 
we shall endow the bundle OV with the measure 

Wohm{d^) = VoLm(dg) ® ay (dm). 

Recall that the point e = (m, g) G OM has been fixed above. 

Proposition/Definition 5. LetYeVe- 

1. The random variable Gfjlff^oo has a smooth density /v(eo ; e) with respect to the measure 
VoLov(c^) on OV. 

2. We have /vK^o ; e) = /v(eo ; e) for any other Y' in Ve- 

So this quantity /v(eo ; e) is independent from Y € Ve/ call it the value at point e of the one- 
particle distribution function of the {V, diffusion started from cq. We shall denote it 
by /(eo ; e); it is defined for e 7^ cq. 

As is clear from its definition, this function takes the same value on points with the same 
HM-projection. Given any point eo G OM, we shall adopt the usual conventions and shall denote 
by 

/+(eo) = /^((?Tio,go)) = {(7(1), g') e OM; 7 future-oriented timelike path, 7(0) = mo, g' G 0^(i) 

the chronological future of gq. This is an open set of OM. It comes from the support theorem 
of Stroock and Varadhan that /(eo ; •) is positive in I+(eo) and null outside the closure of 
/+(eo)a. 

We describe here the proof of proposition/definition [5] without technicalities. 



normal flow lines 




path of a (3, y)-diffusion 



Figure 3: Constructing the one-particle distribution function. 

We use the same idea as in section 12.21 where a family of constant time hyperplanes was 
used to re-parametrize the process. These global objects will be here replaced by local ones: the 
normal variation {Ve}£g(_^ ,j) of the spacelike hypersurface V. Their local definition is illustrated 
in figure 3. Suppose that eo belongs to V_r,; the timelike path {es}s^o will then hit each V^ 



^■'in Minkowski spacetime, this result comes from proposition 8 in the article |19| of I. Bailleul. A similar proof 
can be given in the general framework of Markovian (V^ 3)-difIusions on any OM. 



once, at increasing e. So we can use e as a time parameter in place of s. Using the flow of the 
normal variation, we can consider the re-parametrized path as a hypoelliptic diffusion in V. We 
shall then get the conlcusion from Hormander's theorem on hypoellipticity. 

The second point of proposition/definition [5] is established using the property of the sets 
Um mentionned in the introduction. Indeed, suppose that V' and V not only have the same 
tangent space at m but are equal in a neighbourhood U of m in V. Then, since V is a spacelike 
hypersurface of M contained in Um and since no spacelike path of (M, q) can hit V or V more 
than once, the set of trajectories of the (F, 3)-process hitting OV in OU is the same as the set 
of trajectories hitting OV in OU. So the densities of ejj^, l//^,<oo and ghy'^HyKoo are under Peg 
equal on OU, i.e. /vK^o ; e) = /v(eo ; e) for any e £ U. Shrinking U to {m} formally gives the 
second point of proposition/definition [5l 

< Proof 1. Normal variation of a spacelike hypersurface. Let V G Ve- For m € V 

and e small enough, define 4>e{fn) as the position at time e of the geodesic started from fh 
leaving V orthogonally, in the future direction. Then there exists (as a consequence of the 
local inversion theorem) a positive constant rj and an open set Z// C M such that the map 
(j) : {—rj,'!]) X V — > Z//, (e, m) i-^ 4>£{rn) is a diffeomorphism. It has the following properties, 
where we write for 4>e(y). 

• (j)Q{fh) = fh, 

• decpeim) G H^^(a)M, and 

• decpsirn) is orthogonal to T^^{m)^e- 

The family of spacelike hypersurfaces {Ve}£g(_^ ,y) is called the normal variation of V. 
The open set U has the function e as a time functioiJ^. We shall suppose without loss 
of generality that U is diffeomorhic to an open set of ]R^+'^. The diffeomorphism (p can 
be extended to {—T],r]) x OV — > OU. To that end, given e G {—ViV) transport parallelly 
g G O^M along the path {0i("i)}jg[o e]' '^'^^^^ T'gt-oS ^or the element of T0^(,7j)M obtained 
that way. The map 

(e, (a,g)) G (-r?,r7) x OV ^ (</.,(?n), T/^pg) G OU 

is easily seen to be a diffeomorphism extending 4>. We shall still denote it by cj). 

Notations. Given a point m G V^, we shall denote by w{m) the future unit timelike vector 
orthogonal to T^Ve. We can extend this vector fields on Z// C M to a vector field on OU 
lifting it horizontally; we shall still denote it by w. In addition to this vector field w on OU 
we shall need some more notations. 

• 7 := q{w{e),g^) will be a function of e = (m, g) G OM. 

• The *®'^=-operation will stand for taking the L^(VoLovJ-dual. 

• Last, Hs will denote the hitting time of OV^ C OM. 

Given a point e G OM, we shall denote by I^{e) its timelike past: 

/-(e) = /"((m,g)) = {(7(l),g') G OM; 7 past-oriented timelike path, 7(0) = m, g' G 0^(i)M}. 

The timelike past of a set will be the union of the timelike past of its elements. 

a) We shall suppose first that eo belongs to OU. If eg = (mo, go) does not belong to the 
(closure of the) timelike past of V, then no timelike path started from mo can ever hit V, 



These sets U are those used in the introduction to construct the sets U, 



so the function /(eo ; •) is null in a neighbourhood of OV. As we are interested in what 
happens near OV, we shall make the hypothesis that eg belongs to the timelike past of V. 
We shall suppose, without loss of generality , that eg G V_^; it will be fixed throughout this 
paragraph. 

As the hitting times will be Peo-almost surely finite under the preceding hypothesis, we 
can consider the re-parametrized process {eH^}e£{-T],r])', it has generator L. We shall 
decompose this operator under the form 

Ve = 0,(e), 4^(e) = (tz7/)(e)+L(/o0,)(e) = (t^/)(e) + (I/)(e), (3.12) 

7 

where L is a (smooth) second order differential operator on OV, and where, as a consequence, 
L acts only on OV^. Now, define the OV-valued process {egj^g^^^^^ := {</'J^(e//^)}^g^_^ 

and denote by is its time-dependent generator. The vector fields V and Vi acting only on 
the fibers of the projection OM — > M, it is easily seen that is is a hypoelliptic operator, so 
the random variable has for any e E (— r/, r/) a smooth density with respect to to Volov- 
It follows that en = gq also has a smooth density with respect to Volov- 
6^ To deal with the general case where eg does not belong to OW, we can suppose without 
loss of generality that V is a subset of a spacelike hypersurface V' such that the analysis 
of point a) applies and such that any timelike path hitting OU hits V'_^, before. Then, 
denoting by /i(eo ; e)VoLov'('^) the smooth hitting distribution of OV' by the process e. 
under Peo, we have 

C OV, Pe„ {en eA, H <oo) = (^j h{eo ; e)/(e ; e')VoLov'_^,((ie)^ VoLov(de'); 

from which we conclude that the random variable e//l//<oo has under Pe,) a smooth density 
with respect to Volov, equal to /v(eo ; e) = / /i(eo ; e)/(e; e')VoL(0)v' (de). 

-v' 

2. The formal proof of this point proceeds using a slightly different point of view than 
the heuristic described before the beginning of the proof of proposition/definition El Fix 
e = (m, g) G OM and let r/o > be smaller than the radius of definition of the (Lorentzian) 
exponential map exp„ : TmM — > M, and small enough for the geodesic ball of radius % to 
be included in Um- Given r] < i]q, denote by the hypersurface of M defined as 

Ar, := {exp„(sr) ; |s| < 77, T G (gO)^}. 

For rjo small enough, the hypersurface A^^^ will be spacelike; pick such an 7]q. Denote also by 
the set of points of M of the form ex.p^,{sU) for m' ^ Arj, \s\ < rj^ and U G T^'M. 



2rf 




2r] 



The set B^j has two important properties. We use the notation V for any spacelike hyper- 
surface belonging to Vg. Recall that (Ty stands for the volume element induced by q on 



V. 



„d+2 



ri,0+ 



Cd 



(3.13) 



If we write V^^ for the intersection of V with the chronological past and future of B^j in Um, 
we have 

The constant Cd appearing above is the Euclidean volume of the unit ball of W^. Given now 
any hypersurface V € Ve, < r/ < r/o and a positive integer N, run independent (V,^)- 
diffusions started from Sq. We shall write e^-'^-j for the random position of the i^^ diffusion 
stopped at the random time i/^*^ where it hits OV (provided this time is finite). Associate 
to a given real valued Lipschitz function (p on OUm the random variable 



=1..N 



(3.15) 



The almost sure following limit is a consequence of the strong law of large numbers: 

FNiv) _ 



Af+oo N 



ip{e) /v(eo ; e)VoLov(c?e). 



If we now let H^^ be the hitting time of the set by the i^^ (V^, 3)-diffusion, set 



=l..Af 



Since if is Lipschitz and -B^ has a 'height' of order r/^, we have 



'^Hf <oo 



s; Cif 



for some positive constant C; it follows that 



^'"^ ( — 



N,oo 



Together with equations (|3.13p . (|3.14p . (|3.15p . this equation gives us the existence and the 
value of the limit 



lim ( Vol(S^)3^ lim^^^''^ 



77,0 



N,oo N 



(/9(m,g)/v(eo; (m, g)) VoLm(dg). 



The left hand side being independent of V, the functionnal of the Lipschitz function ip 
defined by the right hand side is also independent of V. The class of Lipschitz functions if 
rich enough to conclude from that fact that the measure /v(eo ; •)V0Lm(-) is independent of 
V € Ve, which implies that /v(eo ; e) itself is independent of V € Ve. > 

To state the next proposition on /(eo ; ■) we shall write V for a spacelike hypersurface of M 
and shall denote by H the hitting time of OV. Given a point e = (m,g) € OV, we shall denote 
by Ci7v(e) the future unit timelike vector orthogonal to T^Y (in accordance with the previous 
notation). This fundamental proposition extends the second point of theorem [2] to the general 
framework adopted in this section. 



Proposition 6. Let eo be a point of OM not belonging to OV. We have 



Eeo [/(eH)lH<oo] = / /(e)(7(g",n7v(e))/(eo; e)VoLov(de) (3.16) 

Joy 

for any bounded function f on OV. 

< Proof - We shall use the notation /v(eo ; •) to denote the (smooth) density of the law of the 
random variable enlfKoo under Pe„, with respect to Volov- Given a point e € OV, we are 
going to prove that 

/v(eo, e) = g(zz7v(e),gO)/(eo; e). (3.17) 

This point e = (m, g) is now fixed. We shall denote by W a hypersurface of Ve; we have 
seen in proposition/definition [5] that /(eo ; e) = /w(eo j s). 

Idea of the proof. The idea of the proof is simple and illustrated in figure 4. Pick a 
positive integer A^; it will be sent to infinity at the end of the proof. Let {Ve}^g^_^^^ be 
the normal variation of V. The positive real rj is chosen in such a way that any timelike 
geodesic started from V_^, of length ^ hits V,,. It implicitly depends on A^; we choose it 
as a decreasing function of N converging to as increases to infinity. We shall write V'_^ 
for the set of points of V_^ from which any future-oriented timelike path hits V. We can 
suppose without loss of generality that the hypersurface W is included in ^ = V^. 

e=(-r?..r;) 







— \ 









V 



Figure 4: Proof of proposition [6l 

Given a point e = (m,g) G OV'_^, we shall write mv(e) for the intersection of the future- 
oriented geodesic 7g started from rh in the direction g'^ and by gv(e) the image at the point 
mv(e) of g by parallel transport along 75. We set 

Xv(e) := (?nv(e),gv(e)). 

The point Xw(e) is defined similarly using W in place of V. Let us denote by the 
hitting time of OV'„^ and set 



Xy:=Xy{&iii <oo e CDV and := ^w(e//' )\h' <oo ^ 



We are going to see that these random points have smooth densities at e which satisfy 
equation (|3.17p : they both depend on N. Equation (|3.17p itself will be obtained as a limit, 
sending A'^ to infinity. 

Proof. Given a (small) open neighbourhood ^ of m in V define ^ C W as the intersection 
of W with the chronological past and future of A in U. Pick A and rj small enough in such a 
way that any timelike path hitting A oi B_ hits V'_^ before. All open sets A used hereafter 
will implicitly be supposed to be included in this fixed A. 

The following lemma is proved noting that the maps Xy : X^^iOA) OA and : 
X^^{OB_) OB_ are well defined smooth diffeomorphisms; as such, the push forwards of any 
smooth measure on OV'„^ by these maps are smooth measures on OA and OB_ respectively. 

Lemma 7. The laws of the random variables XylxyGOA and X^lxweOB under Peo have 
smooth densities with respect to IoaVolov o,nd IobVolow respectively. 
These densities are denoted by /y^^(eo ; •) and f^\eQ ; •) respectiveljl^l. We are going to 
prove that we have 

g(r^7v(e),gO)/(f)(eo; e) = /f ^(eo; e); (3.18) 

we shall then get identity (|3.17p using the following lemma. 
Lemma 8. • f^^\eo;e) /v(eo;e). 

N+oo 

• /vr^(eo; e) /(eo; e). 

iV+oo 

To proceed further and establish identity (|3.18l) . we need to give some definitions. If V and 
r] are chosen small enough, there exists a bundle isomorphism trivializing OU: 

■0 : R^''^ X SOo{l,d) OU. 

We shall denote by (C, g) the point V'"^ (e) = ((m, g)) ; the set Z will be the intersection of 
a small Euclidean ball of Ali+d(M), of center g, with 500(1, d). Set Az := OAnij{R^''^ x Z) . 
The set Bz C W is the defined as the set of points of W of the form expff^{sT^), with 
{m,T) G Az- The collection of all the (m',g') S OBz where m' = exp^(sr'^) and g' is the 
image of T by parallell transport along the geodesic exp^/(- is denoted by Bz. Similarly, 
A'z is defined as the set of points of V of the form exp^(sT'^), with {fh,T) G Bz- The 
collection of all the (m', g') G OV, where m' = exp^(sT'^), {fh,T) G Bz, and g' is the image 
of T by parallell transport along the geodesic expj^j (•r°), is denoted by A'z- The range of s 
in these definitions is restricted in such a way that the geodesies exp.(- T") remain in U. 
To prove identity (I3.18P we start from the inclusions 

{Xy g Az} c {Xw G Bz} c {Xy g A'z} 



to get the inequalities 



■ eo 



[Xy G Az) " Peo(^v € Az) 



I.e. 



^ ^ (ep ; e)lg^(e)VoLow(c?e) ^ Peo (^v € -A'z) ,^ 

" //f ^(eo;e)l^,(e)VoLov(de) " IPeo(^ve^z)' 

We are going to obtain identity (|3.18p taking successively the supremum limit in the above 
inequalities, first as A decreases to {?7i}, and then as Z descreases to {g}. This final step 
rests on the following fact. 



^^Recall 7?, and so Xy and Xw, depend on N . 



Lemma 9. lim lim ") ^ = 1. 

Z\{g} A\{m}Peo(Xv G Az) 

This comes from the fact that the ratio of the Volov_^ -volume of the sets {e € OV'_^ ; Xv(e) G 
^i^} and {e e OV'„^ ; Xv(e) € ^l'^} converges to 1 as ^ \ {m} and Z \ {(/}. Recall ay 
and a"w are the volume element induced by g on V and W respectively. It remains to evalu- 
ate the ratio of the integals in equation (|3.19j) to get the conclusion; to that end we use the 
following fact. 

Lemma 10. • The limit lim ^—-r- exists, and 

A\{m} aY[A) 

• there exists a positive function c{Z) of Z, decreasing to as Z decreases to {g}, and 
such that we have 

(l-c(Z))gK(e),gO) ^^hm^^^^ ^ (l + c(Z))gK(e),g°) 
for Z small enough. 

o To see where this result comes from, write exp^ : T^V V for the exponential map in V 
at point m, and exp^ : T^W W for the exponential map in W at point m. We measure 
volumes in T^V and T^W using the (constant) volume elements Vol^ and Vol^ induced 
by q on T^V and T^W respectively. Writing A = exp^ (^) we have 

"^("^^ ■ 1 (3.20) 



Associate to g G OmMn'0(M^''^ x Z) the set C T^M, image of 2 C T^V by the projection 
map — > T^W parallelly to g". We have on the one hand 



V0Lw(UgBg 



1, 



O'w(-Bz) A\{m} 

where the union is taken over all g G Or„M n X Z) , and on the other hand 



VoC(i?g) 



Together with limit (|3.20|) these two estimates imply lemma [TOl o 

Decomposing Bz into the union of its fibers: Bz =■ B'z ■, we can write the integral 



e) as 



j /if^(eo; e) iBz(e) 

fw\^o; ("i,g))lBffi(g)VoLa(dg) ) lBz{fh)V0Ly^{dm). 



A similar decomposition can be written for J /y^^(eo ; e) 1^^ (e)VoLov('^) using the de- 
composition Az ='■ A^ of Az into fibers: 



meBz 

(TV) 



/v (eo; (m,g))l_4^(g)VoLa(dg) lAz{fn)yoLY{dm). 



Note that Az and Bz have the same fiber at point m, namely = O^M n ipiB}''^ x Z). 
We get as a consequence of lemma [10] the following two inequalities: 

o^/o^M/^fHeo; (m,g))lB™(g)VoLo„M(c^g) //if^(eo; e)li3^(e)VoLow(c^) 

l-c[Z))q[wY{G),^ ) TjyT- ^ hm ^-f^ 

/o™m/v (eo; ("i,g))lBS'(g)VoLo™M(dg) ^\W//v ^^0 ; e)1^2(e)VoLov(c?e) 

and 



r- //if^(eo; e)lB;,(e)VoLow(c;e) ^ , / q /o„m /i^H^o ; (m, g)) lB-(g)VoLo„M(dg) 

A\{rn} j (eo; e)l^^(e)VoLov(de) /o„m /v (^o ; ("^,g))lB-(g)VoLo,„M(dg) 

Taking the supremum limit as Z decreases to {g} and using lemma [10] we obtain 

I fw\^o ; e)lg^(e)VoLow(^) . , . o^ fw\^o'^ ^) 
hm hm - — = q[WY{e),g j-fr-- . 

^\MA\{™}//f )(eo; e)l^,(e)VoLov(de) /f^eo ; e) 

As equation (I3.19P together with lemma [9] tells us that this supremum limit is equal to 1, 
we conclude that 

g(rov(e),g°)/if^(eo; e) = /^^^(eo; e). 
Identity (|3.17ll follows from lemma [8] sending N to infinity. \> 

The property of the one-particle distribution function emphasized in proposition [6] will be 
used to prove the following fundamental theorem. 

Theorem 11. We have L*f{eo ; ■) = in OM\{eo}. 

The approach to relativistic Ornstein-Uhlenbeck process and (general) relativistic diffusions 
developped so far in the work of F. Debbasch and his co-authors relies entirely on a similar 
(manifestly covariant) transport equation, which is given as the fundamental object in their 
approach. This theorem provides a dynamical justification of this approach. 

< Proof - As theorem[TT]is of a local nature, we are going to take for each point e 7^ eo a V G Ve 
and work in the neighbourhood OU of e constructed in the proof of proposition/definition 
[5] using the normal variation of V. We shall use here the same notations as there; notice in 
addition that in an expresion like g(ro(e), g'^) , the vector zz7(e) will be seen as an element of 
TmM. rather than its horizontal lifting. 

The beginning of the proof is exactly the same as in point 1, a) of the proof of proposi- 
tion/definition [5] We repeat it here to ease the reading ; it is quoted between the two stars 

(*)■ 

1) (*) We shall suppose first that eo belongs to OU. If eo = (mcg*^) does not belong to 
the (closure of the) timelike past of U, then no timelike path started from niQ can ever hit 
V, so the function /(eo ] •) is null in a neighbourhood of OV. As we are interested in what 
happens near OV, we shall make the hypothesis that eo belongs to the timelike past of V. 
We shall suppose, without loss of generality , that eo G It will be fixed throughout this 
paragraph. 

As the hitting times -ff^ will be Peo-almost surely finite under the preceding hypothesis, we 
can consider the re-parametrized process {eH^}ee(-ihv)'^ generator ~^L. We shall 

decompose this oeprator under the form 

Ve = 0,(e), :^(e) = (ti7/)(e)+L(/o0,)(e) = (t^/)(e) + (I/)(e), (3.21) 



where L is a second order differential operator on OV, and where, as a consequence, L acts 
only on OVg. Now, define the OV-valued process {ee}^gj._^^^ := {4>7^{'^He)^ ^^(^_^ ^-^ and 

denote by le its time-dependent generator. (*) This operator is seen to be hypoelliptic, so 
the random variable has for any e G (— r/, rfj a smooth density /o(eo ; e, •) with respect to 
VoLov which satisfies the equation 

Ve e {-11, r,), dsp{eo ■,e,-)= ^p(eo ; e, •)• 



e 



stands here for the L^(VoLov)-dual of i^. Let us now denote by VoLjjy the pull-back 
on OV by (ps of the measure Volov^ on OV^, and denote by Gg its density with respect to 

VoLov. Then has a density /i£(eo ;£,■) = — — with respect to Vol^y; it satisfies 
the equation 

dsM^o ■,e,-) + ^M.(eo ;£,•)= ?f'''^/Ie(eo ; •)• (3.22) 



G. 



e 



We have here (.l^^'^g = q^"^^ for any smooth function g. Denote by /i(eo ; £, •) the 
density of e^^ with respect to VoLoVe, and consider /i and G as functions of e and e G V^, 
i.e. consider them as functions defined on the open set U. Then, equation (I3.22P can be 
written 

mfi{eo; ■) + -^fi{eo; ■) = L /i(eo;-)- (3-23) 

The operator L has been introduced in equation (|3.21j) . It is useful at that stage to remark 
that we havS 

— *OVe _ Y* 

as a consequence of the change of variable formula, and since we have a normal variation of 
V. The following lemma is needed to make the final step. 
Lemma 12. We have for any smooth function f 

w*f + wf + ^f = Q. 

o As above, this is consequence of the change of variable formula and the fact that we have 
a normal variation of V. We have, for any smooth function (f with compact support, 

{w*f) (e) (^(e)VoL(de) = J /(e) {w y,){e)VoL{de) = J /(e,e) (9.(/9)(e,e) (e) do (de) de 

= - j (9e/)(e,e)99(e,e) ^^(e) ao{de) de - j {f^p){e,e)^eGe{e) ao{de) de 

= - j (mf + (e) (^(e) VoL(de). 



As a consequence of this lemma we can use the decomposition given in equation (I3.2ip to 
write equation (I3.23P as 

L* (^^^) = 0. 
Proposition [6] enables to conclude that L* f{eo ; ■) = inU. 



^^Recall that VOL is the Liouville measure on OM and that the *-operation is the L^(VOL)-dual operation. 



2) To deal with the general case where eo does not belong to O^, denote by V C V an open 
subset of V such that any timelike path hitting V hits before. Such a manifold V will 
exist provided rj is small enough. For a small enough 5 > 0, the set (/)(_5 5)(V) := {(/)e(m) G 
M ; m G V' and |e| < 5^ will have the property that any timelike path hitting it hits 
before. Set 

Z^':=0</.(_5,5)(V'); 

this is an open set of OM. To prove theorem [TT] on U\ it suffices to remark that for e £U' 
, we have 

/(eo ; e) = E^,, [f{eH ; e)lH<oo] , 



where H is the hitting time of The previous part of the proof applies to each function 
f{ef{ ; •). It follows then from the ab( 

as a mean of L*-harmonic functions. 



f{ef{ ; •). It follows then from the above identity that /(cq ; •) i is an L*-harmonic function 



3.3.2 L-harmonic functions 

We shall see in section 13.41 an important application of theorem [TT] in relation with statistical 
irreversibility. Before turning oursleves to that side, we would like to stress in this section the 
importance that theorem [TT] might have from a geometrical point of view. To that end, we shall 
investigate its meaning in the study of the (0, e.)-diffusion of Franchi and Le Jan. As emphasized 
after definition [3] this (F, 3)-diffusion is the only process of this class determined entirely by the 
geometric background (M, g); this property gives it a special status. Its generator is 

i=l 

We shall call a function on OM satisfying the relation Lf = an L-harmonic function. 
The class of bounded L-harmonic functions and the asymptotic behaviour of the (0, e.)-diffusion 
are two faces of the same object: the boundary at infinity of the manifold (M, g). 



a) Ideal boundaries of manifolds and invariant cr-algebra. Let us illustrate this corre- 
spondence recalling what happens to Brownian motion on some special Riemannian manifolds; 
as the (0, e.)-diffusion, Brownian motion is entirely determined by the geometric environment. 
Suppose (M, q) is a simply connected Cartan-Hadamard manifold: it is a Riemannian manifold, 
diffeomorphic to some M", with curvature bounded by two negative constants. The exponential 
polar co-ordinates (r, 9) G M+ x S'^"^ associated with any point provide global co-ordinates on M. 
These manifolds have the property any sequence of balls {-Bi}ij>o, with constant radius, whose 
centers leave any compact, appear uniformly small when seen from within a compact set: For 
any compact set K and given any e > 0, there exists an index is such that for any i ^ ie and any 
point m £ K, whose accociated system of polar co-ordinates is denoted by (r 9), any point of Bi 
has polar angle 9 contained in a region of of diameter no greater than e. These manifolds 
also enjoy the following property: Given any geodesic {'yt}t^o and any point m', there exists 
a unique geodesic {'yt}t^o started from m' such that the distance between 7^ and 7t remains 
bounded. These properties motivate the introduction of a compactification of M, homeomorphic 
to and where a path converges to a point of tboundary if its polar angle converges in any 

polar system of co-ordinates. Any geodesic converges to some point of the boundarjo- 



See for instance the article |20| of Anderson, or chapter 8 of the book |21| of R. Pinsky. 



From a probabilistic point of view, we can investigate the far end of a manifold looking at 
what happens to Brownian motion {wt}t^o on (M, q) as time goes to oo. To that end, we define 
the invariant a -algebra of the process as being generated by the real- valued functionals F{w,) 
depending only on the asymptotic beahaviour of w,, i.e. satisfying the identity F(^{wt}t^o) = 
F{{ws+t}t^o) for any s ^ 0. 

Back in the above Cartan-Hadamard manifold (M, q), pick a point m G M and denote by 
{rt,9t) the m-polar co-ordinates of wt- It can be proved^^l that {wt}t'^o converges P«,Q-as to some 
random point of 5M, characterized by the fact that 6t —>■ Boo € S"^"^. It can also be shown 
that the invariant cr-algebra is P^„-indistiguishable from the algebra generated by 9^. This fact 
gives a probabilistic meaning to 5M, or, conversely, gives a geometric meaning to the invariant 
(T-algebra of Brownian motion. 

The situation appears to be similar, though subtler, in the Lorentzian framework of Minkowski 
space. Recall the causal boundary C of M^'^ is the ideal boundary of M^'^ characterized by the 
property that two timelike paths {7t}f^o and {7t}t^o converge to the same boundary point iff 
they have the same chronological past: I~ ij) = The following theorem has been proved 

in the articles [H] of I. Bailleul and [23] of I. Bailleul and A. Raugi. It holds for any starting 
point eo of the (0, e.)-diffusion. 

Theorem 13 ([H], [23]). • The W-^'^-part {mj^^o of the {0,e,)- diffusion converges PeQ- 
almost surely to some random point rrioo of C. 

• The a-algebra generated by coincides with the tail a-algebra of {^s}s^o, up to W^g-null 
sets. 

So, we can find back the causal boundary in the probabilistic invariant cr-algebra. This is a 
nice feature that might help clarify geometrically more complicated situations, giving a simple 
probabilistic picture of what happens. J. Franchi has for example undertaken in [18] the study 
of the (0, e.)-diffusion in Godel's spacetime. This space has a trivial causal boundary, reduced 
to one point. Yet, he has been able to prove that the invariant cr-algebra of the process is not 
trivial. This suggested, in return, the definition of a purely geometric boundary. 

b) Poisson and Martin boundaries. The link between geometry and probability illustrated 
above is complemented by the existing link between invariant cr-algebra on the one hand and 
the set of bounded L-harmonic functions on the other hanj^. It is equivalent to determine one 
or the other. The set of bounded L-harmonic functions is called the Poisson boundary of 
(L,M). So, the Poisson boundary of L, the invariant cr-algebra and the geometry at infinity of 
(M, q) may be seen as three faces of a same object. 

Let us give a last picture of the Riemannian/Brownian situation. We shall get a clearer image 
looking at any elliptic smooth second order differential operator Lq on a connected (relatively 
compact) open set D of M". Recall the Martin boundary of {Lo,D) is the collection of 
non-negative Lo-harmonic functions on Martin gave in [24] a methof to construct this 

set and proved that any non-negative Lo-harmonic function can be uniquely represented as the 
barycenter of a finite measure on the set of extreme points of his boundary. This construction 
is now well understood from a probabilistic point of view (see for instance chapter 7 of Dynkin's 

^*See for instance the pioneering article [20] of Anderson, or the article [22] of Y. Kifer. 

^^Consult for instance chapter 8 of the book |21) of R. Pinsky for the Riemannian case, and proposition 8 in 
the article |19| for the hypoelliptic situation appearing in the study of the (0, e.)-diffusion in Minkowski space. 
^°This set contains the Poisson boundary of Lo in D. 



book [25|)- Let us briefly describe it. Denote by G{x,y) the Green kernel of Lq in D, and define 
the function 

K{x, y) = Q^^^^ X eD,y e D\{xo}, 

where xq is a fixed point. Observe that the function K{-,y) is Lo-harmonic on D\{y}, for any 
y € D\{xq}. It follows that we shall construct Lo-harmonic functions on D sending y to the 
boundary of D, provided the limit limyK{-,y) exists. Martin's boundary is made up of all the 
functions obtained that way. A sequence iT/wlw>n of points of D leaving every compact D, and 
such that the function K{-,yn) convergea^^l. is called a fundamental sequence for (Lq, D). In 
short, Martin's theory asserts that the knowledge of fundamental sequences is equivalent to the 
knowledge of the set of non-negative Lo-harmonic functions. 

Pinsky gave in [26] a probabilistic proof of a useful characterisation of fundamental sequences, 
known from potential theoretists beforj^. Lq will denote the L^(LEB)-dual of Lq and {Px}xeD 
the laws of the diffusion {Xt} in D with generator Lq. We shall denote by C its exit time from D, 
and, given a compact subset U of D, we shall denote by Hu the hitting time of U by {Xt}o^t<f . 

Theorem 14 (Pinsky[26]). The sequence {yn}n^o is fundamental for {Lq,D) iff, for any smooth 
compact subset U of D, the sequence of conditional distributions |Pjy„(X/^j^ S • \Hu < C)}„>o 
converges. 

c) A conjecture. Theorem [TT] bringing into play L* and hitting distributions, through /(eo ; ■), 
it is now time to examine it. Note that 

L* = -Ho+\Y.y'- 

4 = 1 

So, L* is the genrator of an OM- valued diffusion analogue to the (0, e,)-process, except that the 
speed = — is past-directed. Call it (0, e ) -diffusion and denote by Peo its law when 

started from eg € OM. It is clear from its construction that the paths of e . started from eo take 
values in the chronological past /~(eo) of eo- Any open set of /~(eo) is visited by the process 
with positive probability, and any spacelike hypersurface V is hit with positive probability. As 
noted in proposition [HI these hitting distributions are determined by the one-particle distribution 

function of the (0, e,)-process. We shall denote it by /(eo; ■), eo G OM. Theorem [TT] can be 
restated as follows. 

Theorem 15. The function / (eo ; •) is L-harmonic on OM\{eo}. 

It is tempting, after reading paragraph b), to renormalize /(eo; •) and try to get possibly 
non-null L-harmonic functions sending the singularity Cq to infinity. This could be made looking 
at 

7(eo; •) 



for some c G OM, or 

(m, g) G 



/(eo; c) 

/ (eo ; (m, g)) 



Uniformly locally on compact subsets of D. 
^^See the notes of chapter 7 in |21j . 



(0, e.)-diffusion 



(O, e .) -diffusion 



Figure 5: (0, e.)-difFusion and (e ., 0)-diffusion. 

As emphasized in proposition [6l the use of the second ratio essentially amounts to look at the 
convergence of the hitting distributions Peol^Hv G • | < oo) of e. on any spacelike (smooth) 
hypersurface V. All this brings us to conjecture the following equivalence. 

Conjecture 16. The following statements are equivalent. 

1. The sequence {e„}„^o is fundamental for L in (M, g). 

2. For any spacelike smooth hypersurface V o/ M the sequence of conditional distributions 

iPe„ (^H^ e-\Hy<oo) 

converges. 

3. For any c G OM, the sequence of L-harmonic function^^ < / ' ^ > converges uni- 

l / (en ; c) J „^o 

formly on compact subsets of /" (c) . 

This fact would explain why the causal boundary of (M, q) is likely to appear in the picture. In 

order for the conditional distributions Pe„(ej:f^ G • | Hy < oo) to converge, the support of each 
of these probabilities has to converge, for any spacelike hypersurface V. This cannot happen 
unless the chronological past /~(e„) = /~((m„,g„)) of e„ converges, i.e. unless the sequence 
{mn}n^o has a limit in the causal boundary of (M, g). Yet, the study of Godel's spacetime 
by J. Pranchi in [18] has made it clear that this geometric boundary might not be appropriate 
to describe the Poisson or the Martin boundary in some situations. Note, yet, that the above 
analysis using the one-particle distribution function does not apply in this non-strongly causal 
spacetime; no good definition of one-particle distribution function is available at the moment in 
such a framework. 

Last, we should oppose the difficulty of this problem on the large scale structure of (M, q) to 
the previously mentionned fact that the local geometry of spacetime can be recovered looking at 
the pathwise behaviour of the (0, e.)-process. Complications come from infinity... We shall come 
back to the above conjecture in a near future. 

^^Given a compact set K and a sequence {e„}„^o of points of OM leaving every compact, the function / (e„ ; ■) 
is well defined on K for n large enough. 



3.4 H-theorem 



We give in this last section a proof of the analogue of the H-theorem first proved in [27] for the 
R.O.U.P. in Minkowski space, as defined there through a 'Kolmogorov equation'. It has been 
then extended in [28] to the R.O.U.P. (as defined in [16]) in any Lorentzian manifold, and finally 
in [5] to a larger class of 'diffusions' in Minkowski space. We deal here with the general case 
of (y,3)-diffusions in any Lorentzian manifold. We refer to the articles [27], [28] and [12] for 
physical motivations. 

Let U C OM be a relatively compact open set and / and g be two positive smooth functions 
on U satisfying the relations L*f = L*g = 0([fl)- We shall denote by F a continuous unit vector 
field on U; define the function p : OU — > (0,+oo), (m,g) i— > q{Ym,g^)- 

We shall make the following assumptions on / and g; they are sufficient to ensure the existence 
of the integrals below, and to differentiate them. 



• In ^ is bounded. 



• There exists positive constants C and e such that / and its first and second derivatives are 
uniformly bounded by Ce~^ '("i->g) U. 

Define now on U the vector field 



X{m) = - [ g° /(m, g) In 4^ VOL„(dg) 
JomM 9{m, g) 



The main result of this section is the following theorem; no assumption on the geometry of space 
or on the data V,^ is needed. 

Theorem 17 (H-theorem). We have diyX ^ for any two L* -harmonic functions f and g, 
and X defined as above. 

We shall begin the proof of this theorem proving the following lemmsl^. 

Lemma 18. Given any (good) smooth function h on OM set 



Then, 



X{m)= g°/i(m,g)VoL„(dg). 
(divX)(m)=/ {HQh){m,g)\o^m{dg). 



< Proof - Given a path 7 in M and two time s,t we shall denote by T]^^ : T^.M T^,M 
the parallel transport operation along the path {^r}re[s,t]- It is an isometry between the two 
tangent spaces. We shall denote by V the Levi-Civita connection on (M, g). Recall that we 
have 

T^7 Y — Y 

y X = hm 

s^O S 

Recall also that the divergence of X is the (Lorentzian) trace of the map V.X. It means 
that given any choice of orthonormal frame g of T^M, the sum 



d 

i=0 



(z(r,r)(/(v^,x,r) 



^*These functions could for example be of the form /(eg ; ■) and /(eg ; ■) if the strong causality assumption on 
M., q) is satisfied. 

^'Compare with the Appendix to the article [28] of F. Debbasch and M. Rigotti. 



is independent of g € O^M; this is, by definition, (divX)(m). Last, recall that the vector 
field Hq is defined as the generator of the lift to OM of the geodesic flow on HM. Its 
dynamics {(ms,gs)} is determined by the condition = and the fact that gs is 

parallelly transported along the path {ms}. 

Choose now a frame g G O^M and a path 7* such that 7*(0) = m and 7*(0) = g*. Since 
parallel transport is an isometry, we can write 

X^. = [ gO h{Ys,g) VOL^. (dg) = / r,togO h{Ys, TtoS) VOL„(dg), 

TtsX^l = [ g° h{jl, TtoS) VoL^(dg) 



SO we have 



and 



s 

Send s to and sum over i to get the result: 

d 



(divX)(m) = V(7(r,r)<7(Vg,X,r) = / (//o/i)(m,g) VoL„((ig). 

With this lemma in hand we can prove theorem dll Recall that L = Ho + V+ \ViB'Wj and 
L* = -Hq + V* + \ViBWj. 

< Proof - First, use lemma [TSl to write 

-divX = j (/In^ VoL(dg) = j (Hof) In^ VoL(dg) + j (Fq/ - ^Hog) VoL((ig). 

Use then the relations L* f = L*g = to get 
- divX = j (V*f + lVi{B'Wj{Xf))) (in f + 1) VoL(dg) - / f (V*g + ^V,{B^^V,iXg))) VoL(dg) 

= / {{V*f){lnl + 1) -^V*g) VoL(dg) + 1 / ((inf + l)y,(i?^^Y,(A/)) - ^V,{BW,{Xg))) VoL((ig). 

Integrating by parts and using the relation y(ln^) = the first integral is seen to be 

equal to 



Recall V* = —Vi. Use integration by parts in the second integral and the relation Vj{Xf) 
V,iXg)j + XgV,[{),to get 

-divX = ^ j y,(ln^)s'^y, (A/)VoL(dg) + ^ j y,(0i?^^yj(A5)VoL(dg) 

= (^^V,(^)B^^V,iXf) -V,[^)B^^V,{Xg)^ VoL(dg) 



2 



- / {H^B'^ {y,{Xg) + ^V,{t^^-V.{^)B^^V,iXg)'^ VoL(dg) 



(3.24) 



We get the conlustion from the non-negativeness of the matrix B = [A ^ ^ ■ 

• Attention should be paid to the range of apphcation of theorem [T71 It seems tempting, 
indeed, in a strongly causal spacetime, to apply it to functions of the form /(eo ; ■) and /{e^ ; •). 
However, today's state of art is far from being sufficient to provide estimates on these functions 
good enough to ensure that hypothese like those made at the beginning of the section holj^. 
This is a difficult topic where the non-ellipticity of the operator L complicates everything. As a 
first step towards such results, it would be interesting to determine small time estimates of its 
heat kernel; known results are unsufRcient to answer this question. 

• Theorem [17] proves that the flow of the vector field X is volume increasing. It is not clear 
how one should interpret this result from a physical point of view when M is different from M^'"^. 
In this special case, choose a rest frame g G S0{1, d) and denote by t its associated time. Then, 
the integral of X over any hyperplane of constant time is an increasing function of t (provided X 
is equal to at space infinity). This fact justifies that we should call theorem [17] an if-theorem 
in that case. Things are less clear in any Lorentzian manifold, where time does not exist globally. 

Things are even less satisfying from an information theoretic point of view. Recall that the 
relative entropy of a probability P with respect to another probability Q is infinite if P is not 
absolutely continuous with respect to Q, and equal to 



In — — 



if P is absolutely continuous with respect to Q. We write Ep for the expectation operator 
associated with P. Relative entropy is always non-negative, as is clear from the inequality 
a In I ^ a — 6. 

Suppose for clarity that P and Q are probabilities on [0, l] (p^ ). Let Xi, X2, ... be i.i.d. random 
variables, with common law P or Q. Then, given any numbers xi, in [0, 1], we have 



P'^°° will stand for the product measure P (g) P (g) • • • on [0, 1]^. Taking now the Xj's to be i.i.d. 
random variables with common law P, it follows from the law of large numbers that we have 
P®~-almost surely 

Q®n(Xi =Xi,...,A:„ =X„) n,oo 

in a sense that should be made more precise. The above estimate roughly means that the support 

of the probability P®" in [0, 1]" has Q®'^-measure of order e "'^^t^" do] ^ when n is large. 

One owes to the statistician Charles Stein a rephrasing of this fact in terms of tests, which 
should be clear from the above descriptiorEl. 



Lemma 19 (Stein). Let he i.i.d. [0, l]-wa/«ed random variables with common law 

P. Consider the hypotheses "Hq : P = P and "Hi : P = Q and suppose we want to test 
hypothesis Hq against Hi. The quality of a decision region An C [0, 1]" is measured by the errors 



■^^The boundedness hypothesis on |- is even most likely to be untrue for such functions. 

^'^This is not a serious restriction as any probability on a Borel space is isomorphic to a probability measure on 
[0, 1]. This class of spaces is large enough to encompass most of the useful situations. See the Appendix of the 
book |29) of Dynkin and Yushkevich. 

Consult for example section 11.7 of the book of Cover and Thomas for a proof of this lemma. 



P®"(^^) and Q®"(A„). Given e > 0, set = inf{Q®"(A) ; C [0,1]", P®"(^^) < e}. 
Then we have 

limlm = -H(F; Q). 

£,0 noo n 

To understand this lemma, imagine you want to test the hypothesis "Hq : P = P", with a 
given (very) small bound on the two errors. Then, the smaller H(F ; Q) will be, the bigger n 
will have to be in order to design a test achieving the requirements on errors. 

If now P and Q depend on some 'time' s and H{¥s ; Qs) decreases, then you will need more 
and more data to achieve the test, as time passes. It could be said of a situation where H{Fs ; Qs) 
decreases to that the process X2(s)...) forgets its law as time increases, as it is more 

and more difficult to distiguish if it has common distribution Pg or Qg. 

From that point of view, a satisfying -ff-theorem for the (V, 3)-processes would be completely 
different from theorem [iTl Given a path 7 in M and a (proper) time s (of 7), define 

. ^_ 9(7^,g°)/(eo; (7s,g)) 



/o^,M9(7s,g°)/(eo; (7s, g)) VoL^,(dg)" 

Let V be any (small) spacelike hypersurface such that 7^ E V and (7s)''" = T-y^V. The probability 
h2{eQ] (7s,g))VoL^^(dg) is the conditional law of the random variable ghy'^HyKoo, given that 
n^HY = 7s; we have seen in section [3.3.11 that this conditional probability does not depend on 
V but only on 7<j and 7^. Given two initial conditions eo,eQ of the (F, 3)-diffusion, define the 
time-dependent relative entropy associated to the path 7 as 

^:o.e(,(^) ■■= I /^Keo; g)ln j;^|5^V0H(dg). 
° JO^.M IT'S (,eo , gj 

A satisfying //-theorem would take the form of the following conjecture. 

Conjecture 20. The MM-valued {V,}) -diffusion process (7., 7') = {(ms,g°)}^>Q e HM 
almost surely forgets its law as time increases: the relative entropy Hj^^ ^/ (s) decreases Fe^-almost 
surely to 0, for any eo,eQ G OM. 

It would also be interesting to see if the following holds. 

Conjecture 21. A freely falling observer has more and more difficulties in distinguishing Peg 
from Fpi . 

We shall adress these questions in a near future. As a last comment, let us notice that theorem 
[17] can be given an information theoretic fiavour. Consider indeed that each open set of spacetime 
initially has a quantity of "information" equal to its volune, and that this "information" travels 
with the fiow of the vector field X. Then, theorem [T7] means that the quantity of information 
that can be found in a fixed open set decreases as the fiow-time increases. Yet, this interpretation 
is far from being as clear as the above two conjectures. 



4 Comments 

It is now time to forget the details of the proofs and summarize the main ideas and results 
exposed above. 

A general class of relativistic diffusions was first presented in the article [1]. Although the 
authors only consider dynamics in Minkowski spacetime, their class of processes is essentially the 



same as the above class of (y, 3)-diffusions. It is characterized by the existence at each time of 
a rest frame with the property that the moving object has, in addition to a deterministic acce- 
learation, a Brownian acceleration in any spacelike direction of the rest frame, when computed 
using the time of the rest frame. Yet, the authors' analysis of the situation rests entirely on a 
transport equation; an approach similar in spirit to the semi-group analysis of Markov processes, 
as opposed to the pathwise study of the process. We propose in this article a simple and direct 
construction of relativistic diffusions on any Lorentzian manifold as flows of stochastic differen- 
tial equations. This construction necessitates to build the diffusions in the orthonormal frame 
bundle of (M, g), as was done by Malliavin or Elworthy for Brownian motion in a Riemannian 
manifold, and by Pranchi and Le Jan in the Lorentzian framework. This change of framework 
is worth being made. Not only are we able to recover directly many of the results established 
so far, but this pathwise approach presents several other advantages over the analytical method 
used up to now. 

First, it provides a direct (co-ordinate free) description of the dynamics in OM which is 
given as the fundamental mathematical object of the model. Simple hypotheses can be given 
(section 13. 2p to construct a diffusion in the more familiar phase space HM from the diffusion 
on OM. An interesting outcome of this approach is the clear new definition of the one-particle 
distribution function that can be given using the pathwise behaviour of the (F, 3)-process (propo- 
sition/definition [5]). The fundamental equation it satisfies (theorem [TTI) provides a dynamical 
justification of the approach used up to now, and sheds some light on the study of the Poisson 
and Martin boundaries of the (0, e )-process (section [3.3.2p . Last, but not least, the formalism 
of vector fields enables us to give in section [3^ a concise and clear proof of a general iJ-theorem. 

Although these results are encouraging, it would be desirable to discuss the adequacy of the 
models provided by (y,3)-diffusions to situations of physical interest. To paraphrase what was 
written in the introduction of the seminal article |^, the models provided by (V. 3)-diffusionJ^ 
should not be considered as accurate models of motion of a "colloidal particle immersed in a real 
(relativistic) medium". Rather, they should be considered as toy models designed to provide a 
framework for the study of the main characteristics of the diffusion phenomenon. In this direction, 
it would certainly be useful to develop an approach to the relativistic Boltzmann equation in 
terms of hydrodynamic limit of a system of interacting particle^. Propagation of chaos results 
could justify the use of (y,3)-diffusions as models of diffusion dynamics. Other dynamics, as the 
one introduced by L. Markus in the article [32], might happen to be of physical relevance. Note 
also the interest that the possibility to define (F, 3)-diffusions in non-isotropic media might have. 

Nevertheless, one can consider as one of the merits of our approach the fact that it provides 
new questions. A few of them have been written under the form of conjectures in sections [3.3.21 
and 13.41 we would like to put forwards two other problems concerning the (0, e.)-process, as we 
think this is a fundamental object. 

Lifetime. The question of explosion of a general (F, 3)-diffusion may appear irrelevant from 
a physical point of view, after reading the above comments. Yet, the study of this problem 
for the 'geometric' (0, e,)- diffusion might happen to be extremely fruitful in its possible links 
with the existence of singularities of the spacetime itself. Indeed, all the studies made so far, 
in Minkowski, Robertson- Walker, Schwarzschild and Godel spacetime^ tend to reinforce the 

''''By the R.O.U.P. in this article. 
Consult the article |12) for a discussion of Boltzmann equation in a relativistic framework. The article |31| of 
Andersson and Comer is also a valuable source of information on relativistic fluid dynamics from a macorscopic 
point of view. 

^^In [IS], [23], [3] and [18]. 



feeling that the M-part of the (0, e,)-diffusion eventually behaves like a lightlike geodesic!^. So 
its seems natural to ask the following question. 



Open Problem 22. Is null geodesic incompleteness equivalent to explosion of the {0,e^) -process 
with positive probability? 

This link between geometry and probability would provide a new approach to the existence 
of singularities on Lorentzian manifolds. It would be interesting for instance to see wether the 
hypotheses of Penrose's theorenJ^ are relevant from a probabilistic point of view or not. One 
of its potential benefits is that the explosion problem has an analytical counterpart which is a 
linear problem. Explosion is equivalent to any of the following two conditionj^. 

1. Let A > 0. There exists a non-null bounded smooth function / such that {L — A)/ = 0. 

2. Let T > 0. There exists a non-null solution to the Dirichlet problem dtf = Lf, on 
[0, T] X M, with initial condition 0. 

The use of the one-particle distribution function /^(eo ; •) of the (O, e ,)-process killed at constant 
rate A will certainly help to see if condition 1 holds. To begin with, it would be interesting to 
find an example of a geodesically timelike complete Lorentzian manifold whose (0, e,)-diffusion 
explodes. No such manifold has been found yet. 



A probabilistic interpretation of Einstein tensor? In so far as the local geometry of 
spacetime can be recovered from the pathwise behaviour of the (0, e.)-proces^, it is tempting 
to ask if one can ultimately give a probabilistic interpretation of Einstein tensor determining 
matter in terms of (0, e.)- diffusion. This question brings us far from the present day knowledge... 
We hope it will have some day a positive answer. 

Aknowledgements. I would like to thank Jacques Pranchi for his numerous comments on an 
early version of the manuscript; they led to a clearer exposition of the results exposed here. 
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